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XXVII.BAND SECHSTES (SCHLUSS-)HEFT 1960 


Die Bewegungsgesetze des Raumkompasses 
Von W. Bauersfeld* 


Richard Grammel zum siebenzigsten Geburtstag gewidmet 


1. Einleitung. Unter dem RaumkompaB ist eine von Max Schuler angegebene Variante des 


Anschiitzschen Kreiselkompasses zu verstehen, die theoretisch gegen Beschleunigungen véllig 


unempfindlich ist. 
Bei der nachfolgenden Behandlung dieses Kreiselkompasses werden ausgiebig Vektoren benutzt. 


_ Es zeigt sich dabei, daB man aus einfachen Ansatzen durch reine Vektorrechnung auf ganz kurzem 


Wege einen iiberraschenden Einblick in alle Bewegungsgesetze des Raumkompasses bekommt. 
Dabei werden nur die kleinen Drallwerte, die infolge langsamer Drehungen der Kreiselgehause 


entstehen, gegeniiber dem Drall der Kreisel selbst vernachlissigt. Ferner werden die Drallachsen 
_ der Kreisel als mit ihren geometrischen Achsen zusammenfallend behandelt, wie es bei technischen 


Rechnungen mit schnellen Kreiseln zulassig und iblich ist. 


2. Die Grundgleichungen. Die Konstruktion des Raumkompasses ist in Abb. 1 schematisch dar- 


_ gestellt. In einer vollig geschlossenen Schwimmkugel, mit der ein Dreibein q,, dy, A, fest verbunden 


oN 


soll durch den Ortsvektor t, die ihres Schwer- 
_ punktes durch 1, dargestellt werden. Zwischen 


verschoben sein. 


der Mittelpunkt der Erde sein, die bei dieser 


wird, sind zwei sehr schnelle Kreisel in zu a, a, 
parallelen Prazessionsachsen drehbar gelagert und | 
durch Zahnrader so miteinander verbunden, da 
ihre Drallachsen }, und ), mit der Richtung 
von Q, gleiche Winkel einschlieBen. 

Die Lage des Mittelpunktes der Schwimmkugel 


beiden soll die Beziehung bestehen 
t= 1, + 0;¢. (1) 


Der Schwerpunkt soll also in Richtung von — q, 
um die Strecke a gegeniiber dem Kugelmittelpunkt 


Der Nullpunkt der Ortsvektoren 1 und 1, soll 


Rechnung einfach als eine Kugel mit dem Ra- Abb. 1. Schematischer athe des Kreiselkompasses. 
dius r behandelt wird. Dabei sollen aber t und 1, 

nicht als mit dem Erdkérper fest verbunden angesehen werden, sondern sie sollen auf einen dre- 
hungsfreien Raum bezogen werden, in dem die Richtungen nach den Fixsternen als zeitlich unver- 


- anderlich zu betrachten sind. 


Zuerst soll der Schwerpunktsatz fiir das Gesamtsystem aufgestellt werden, dessen Masse mit m 
bezeichnet wird. Als auBere Krafte kommen nur in Betracht das Gewicht — (t/r) mg und eine 
Stiitzkraft @, die von dem allseitigen Fliissigkeitsdruck auf die Schwimmkugel herriihrt, und die 
daher genau durch den Kugelmittelpunkt geht. Denn die Reibungskrafte, die bei den langsamen 


z Drehungen der Kugel gegeniiber der Schwimmfliissigkeit entstehen, kénnen vernachlassigt werden. 


Wor 


Damit liefert der Schwerpunktsatz die Gleichung 


a ib 
= Q—_mg 
1b. m R F 5 


‘oder, in Verbindung mit (1), 


& . Yt 
& =tm—a,am+-—-mg. (2) 


* Gestorben am 28. Oktober 1959. 
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Weiter soll zunichst der Flachensatz auf das Gesamtsystem angewandt werden, mit dem Schwer- 
punkt 1, als Bezugspunkt. In diesem Fall tritt als auBeres Kraftemoment nur das von der Stiitz- 
kraft @ herriihrende auf, mit dem Hebelarm a, a. Der Gesamtdrall setzt sich aus den beiden 
Kreiseldrall-Werten 0,B und },B zusammen. Damit wird 

(d, +d.) B= aa, x & 
oder mit (2) 
(+8) B= man, x (tia ttt). (3) 


SchlieBlich soll noch der Flachensatz auf die beiden Kreisel einzeln angewandt werden. Die ” 


Cehiuse dieser Kreisel sind drehbar in zu a3 parallelen Prazessionsachsen gelagert. Sieht man diese 
Lagerungen als reibungsfrei an, so kénnen die Mo- 
mente, die von den Lagerkraften herriihren, keine in 
die Richtung a, fallende Komponente aufweisen. Sie 
kénnen also durch die Ausdriicke a, x %, und a, xX 5 
dargestellt werden. 

Weitere Momente entstehen durch den Zahndruck 
zwischen den beiden Zahnradern, die, wie schon ge- 
sagt wurde, die beiden Kreiselgehause so mitein- 
ander verbinden, daf ihre Ausschlagwinkel ¢ (Abb. 1) 
sich stets gleichen. Man erkennt leicht, daB die Vek- 
toren, die diese Zahndruckmomente darstellen, ein- 
ander gleich sein miissen, und daf sie auch die gleiche 
Richtung (a;) aufweisen miissen. Sie kénnen also in 
beiden Fallen durch a, Z dargestellt werden. 

SchlieBlich soll noch vorausgesetzt werden, da 

Abb. 2. Schema der Federanordnung. die beiden Gehduse, die die Kreisel umschlieBen, durch 

Federn (Abb. 2) so miteinander verbunden sind, daB 

weitere Momente entstehen, die vom Ausschlagwinkel ¢ abhingig sind. Diese beiden Momente 

sollen sich gegenseitig aufheben, kénnen also durch — a3 F und + a; F dargestellt werden. Damit 
entsteht das Gleichungspaar 


Pee ean cee ye 


db, B= a, xX, +03 F + 0,2. 


Hierzu kommen noch die einfachen geometrischen Beziehungen fiir die Einheitsvektoren }, 


D, — Dp = 20, cose. 

3. Die Hauptgleichungen fiir die Bewegungsgesetze. Aus den Gleichungen (4) und (5) folgt sofort 
(d, + d,) B=2a, Bsine + 20, Be cose =a x (4, + &) +2052, 
(, — b,) B = 24, Bcose —2 a, Be sine =a, x (CQ, —&)—20,F. 


Diese beiden Gleichungen werden jetzt skalar mit @; multipliziert. Da nach (3) ay+(D, + d,) B=0 
werden mu, so entsteht 


0=2,-a; Bsne = 2Z 
und 
20:0, Boose =—2F. 


Hiernach wird der Zahndruck Z zu null, und, da wegen. (lg: 03 = 0 und a,- a, = 0 


(2° 3 = —,- Ay und i, * Os = — 0+ Oe 
1st 


A: 3 = 0 und i) = 
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Die Giltigkeit der ersten dieser Gleichungen ist an die Bedingung ¢ > 0, die der zweiten an 
€< 7/2 gekniipft. Nun ist aber, ganz allgemein, 


Gg = My Ay + Os + Ay M+ Og + A3 Ag+ dg» 


und da wegen (3-3; = 1 auch aj: a3 = 0 ist, so gelangt man sofort zu der ersten Hauptgleichung 
dieses Kreiselsystems 


. : F 
Daas peli Bae Oe axe (6) 


Sie besagt, daB die a,-Achse sich stets in die Richtung von Qz einstellt. 
Die zweite gewinnt man nun auch sehr schnell. Nach (6) ist 


PS SU aa erer 
und damit nach (5) 


: - Bsi 
D, 4 Dy = 2:0, sin € =a, Se eee 


Die Einsetzung dieses Ausdrucks in (3) liefert 


é : bes Asta ae nd) BA sineare 
(Bs + By) B= a5 x tp —E* + 05 X Ga (“| 


=mady x (ii a +14) 


oder, mit Einfiihrung der HilfsgréBe 


B sin 2) ¢ 


ee eta ae? @) 


x fe + Sa, + «) tnd} 0. (8) 
Das ist die zweite Hauptgleichung. 


4. Die Schulersche Lésung. Die Gleichung (8) wird wesentlich einfacher, wenn die HilfsgréBe 0 
zu einer Konstanten wird. Dazu ist nur nétig, wie (7) zeigt, da} die Federkonstruktion zwischen 
den beiden Kreiselgehdusen so ausgefiihrt wird, dai das Moment F' aus den Federkraften pro- 
portional sin 2 ¢ wird. Die in Abb. 2 schematisch eingezeichnete Federkonstruktion erfillt diese 
Forderung. Die beiden Federn sind véllig gleich und in der Stellung ¢ = 0 auch gleich gespannt. 

Wenn man nun weiterhin die in (7) enthaltenen Konstruktionsdaten B, M, A und F = x sin 2 € 
so wahlt, daB die HilfsgréBe 0, die die Dimension einer Lange hat, gleich dem Erdradius r wird, 
so laBt die Hauptgleichung (8), die in diesem Fall die Form annimmt? 


a, x(E+1 4 ir) =0, (9) 


eine tiberraschend einfache Lésung zu. Sie wird namlich mit dem Ansatz 


a (9a) 
identisch erfillt. Das bedeutet: Wenn die Achse a, zu Anfang genau die Richtung von rt hat, also 
die durch den Erdmittelpunkt bestimmte Richtung, so behalt sie diese bei beliebigen horizontalen 
Geschwindigkeiten t und Beschleunigungen { starr bei. Zu solchen horizontalen Geschwindig- 
keiten ¢ gehért auch die durch die Erddrehung bedingte. 

Diese iiberraschende Lésung des KreiselkompaBproblems stammt von Max Schuler. Der Weg, 
den er hierfiir beschreiten muBte, war sehr viel miihsamer als der vorstehend benutzte, der relativ 
leicht zu finden war, weil das Ziel festlag. Die Bezeichnung ,,Raumkompaf« ist fiir dieses System 
von der Firma Anschiitz in Kiel eingefiihrt worden, bei der Schuler damals als wissenschaftlicher 
Mitarbeiter tatig war. 


1 Der Schwerpunktsabstand a, der nur die GroBenordnung von Zentimetern hat, fallt natiirlich in Gl. (8) 
gegentiber der Lange g = r ganz heraus. 
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Mit dem Ansatz (9a) scheint zunachst die Schuler-Lésung auf horizontale Geschwindigkeiten 
beschrankt zu sein, mit ¢ senkrecht zu r. Denn nur fiir solche bleibt 1/r ein Einheitsvektor, da r 
als ein fester Wert fiir den Erdradius eingefiihrt worden war. Diese Beschrankung besteht aber 


praktisch nicht, wie spater gezeigt werden soll. 


5. Lotanzeige und Nordweisung. Der RaumkompafB stellt theoretisch ein ideales Gerat fiir die 
Lotanzeige auf einem schwankenden Schiff dar, weil seine a,-Achse standig nach dem Erdmittel- 
punkt ausgerichtet bleibt. Die Anzeige der Nordrichtung geht aus der ersten Hauptgleichung 
hervor. Diese Gleichung (6) besagt, daB die Richtung der Kugelachse a, stets mit der von (lg Zu- 
sammenfillt. Fiihrt man den Einheitsvektor ) fiir die Richtung der Erdachse ein, so wird @, fir 
den auf der Erde fest aufgestellten Raumkompa8 dargestellt durch 


ls = p X Ag Wo (10) 
mit W) = 7.29 - 10 sec fiir die Umdrehungsgeschwindigkeit der Erde. Daher wird 
3 = Ay Wp COS Po » (10 a) 


wenn der Winkel zwischen )) und a3 = t/r mit (7/2 — qo) bezeichnet wird. Damit kennzeichnet 
die geometrische Polhéhe, die wegen der Zentrifugalbeschleunigung der Erde etwas kleiner ist als 
die physikalische Polhéhe y, die wir auf der Erde messen. Der Unterschied zwischen beiden be- 
tragt, wie man leicht feststellen kann, nur 

ie 


2g ow; sin 2m = 5.94sin2q in Winkelminuten. 


Da die Richtung der Achse q, identisch mit der von ) X az ist, fallt sie genau nach Osten, und 
daher die von a, genau nach Norden, mit dem kleinen eben angegebenen Winkel als Erhebung oder 
,Blevation tiber den Horizont. 

Auf dem fahrenden Schiff erfahrt t eine von der Kursrichtung und der Schiffsgeschwindigkeit 
abhangige Richtungsanderung, der natiirlich auch Q, unterliegt, und die eine als ,,Fahrtfehler* 
bekannte Fehlweisung von qa, und a, hervorruft. 

Aus der ersten Hauptgleichung (6) gewinnt man weiter den Winkel ¢, der sich im RaumkompaB 
einstellt. Nach (7) wird ja fir 9 = r das Federmoment 

2 


B 5 
(hs sin2e. (11) 


ar 


Steht nun der Raumkompaf fest auf der umlaufenden Erde, so wird nach den Gleichungen (10a) 
und (6) 


F 
ioe 10'S Beoaee 
und daher in Verbindung mit (11) 
. a mar Mo COS Po 
sin € = > B ; (11a) 


Diese Gleichung ist fiir die Konstruktion des Raumkompasses auSerst wichtig. Denn ¢ darf ja 
den Wert z/2 nicht erreichen. Fiihrt man noch den Winkel ¢, ein, der sich nach (11a) am Aquator 
(% = 0) einstellt, so gilt 

2 Bsiné) = mara,. (11b) 


Man ersieht hieraus, wie man durch passende Wahl des Schwerpunktsabstandes a den Winkel ¢ 
konstruktiv festlegen kann. Mit diesem laBt sich die erste Hauptgleichung (6) auf die pinfache 
Form bringen 
raat Oren 
As = 0, Oo = (12) 

6. Die Eigenschwingungen des Systems. Die Gleichun ie fii 

‘ g (9), die fir den RaumkompaB di oi 

Hauptgleichung darstellt, gibt nun auch AufschluB iiber die Eigenschwingungen des S stan ante 
dem Ansatz (9a) wurde ja nur eine singulare Lésung gewonnen. 

Die médglichen Bewegungen, die Gleichung (9) darstellt, lassen sich am einfachsten tibersehen 
wenn man diese Gleichung mit der allgemeinen Bewegungsgleichung des F adenpendels vergleicht, 
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Um diese aufzustellen, soll die Pendellange mit / bezeichnet werden, die Masse des kleinen Pendel- 
kérpers mit m. Ferner wird zur Bestimmung der Lage des Aufhangepunktes der Ortsvektor r 
und fiir den Schwerpunkt des Pendelkérpers der Ortsvektor 1, eingefiihrt, wobei t und t. in genau 
der gleichen Weise wie bisher auf den Erdmittelpunkt bezogen werden. Die Richtung des Fadens 
soll durch den Einheitsvektor a, bestimmt werden (Abb. 3). Daher ist 


tr, =t—agl. (13) 


- Die Krafte, die am Pendelkérper angreifen, sind die Gewichtskraft — (t/r) mg und die Zugkraft 
des Fadens, die mit § bezeichnet werden soll. Nach dem Schwerpunktsatz wird daher 


t,m = R — = mg. 
In Verbindung mit (13) ergibt sich daraus 
RK = m(i— Gig +1 a 
oder, da § die gleiche Richtung wie a, hat, 


ay x (E+ 1 £— fil) = 0. (14) 


Diese Gleichung stimmt mit der Gleichung (9) des Raumkompasses viéllig itberein mit Se Oe aa 
es einfachen 


der einzigen Abweichung, daB in dieser als Faktor von (jz der Erdradius r anstelle Fadenpendels. 
der kurzen Pendellange | auftritt. 

Daraus ist zu schlieBen, daB der RaumkompaB um seinen Stiitzpunkt (r) mit seinem Schwer- 
punkt (r,) ganz gleichartige Bewegungen auszufiihren vermag wie das einfache Fadenpendel, 
namlich ebene Schwingungen, Kreisbewegungen und ellipsenartige Schwingungen. Aber diese 
Bewegungen verlaufen auBerordentlich viel langsamer als beim Fadenpendel. Wahrend das Faden- 
pendel bei kleinen Ausschlagen die aus (14) hervorgehende bekannte Schwingungsdauer 


ee) n/t (14a) 

& 
von der Gré®enordnung von nur wenigen Sekunden hat, wird sie analog beim RaumkompaB zu 
{As yz = 84,4 Minuten. ; (14b) 


Die Bewegungen des Raumkompasses verlaufen hiernach im Verhaltnis 1:2500 langsamer als die 
eines Fadenpendels von 1 m Linge. 

Aus dieser Analogie erkennt man, da der Raumkompaf auch gegen vertikale Beschleunigungen 
unempfindlich ist. Denn ein ruhendes Fadenpendel kann man nicht durch Auf- und Abbewegen 
seines Aufhangepunktes zum Schwingen bringen. Nur dann, wenn es schon anfangs hin- und her- 
pendelt, kann man die Schwingungen steigern, falls man bei der Auf- und Abbewegung einen be- 
stimmten Rhythmus innehalt. Da der Raumkompaf} aber durch horizontale Beschleunigungen 
nicht aus der Ruhelage herausgebracht werden kann, miissen auch vertikale Beschleunigungen 
wirkungslos bleiben. 

Nun entsteht weiter die Frage: Wie verhalten sich die Achsen qa, und a, der Kreiselkugel sowie 
der Winkel ¢, falls die a,-Achse wie ein Pendel schwingt ? Dariiber gibt die erste Hauptgleichung (6) 
Auskunft. Diese Gleichung war bereits mit dem Federgesetz des Raumkompasses auf die Form (12) 
gebracht worden. Sie besagt zunachst, daB die Kugelachse a, in jedem Zeitpunkt die gleiche 
Richtung wie Q; hat. Wenn also z. B. bei der Aufstellung des Raumkompasses am Nordpol diesem 
eine solche Schwingung aufgepragt wiirde, daB die Pfeilspitze von a, eine Kreisbahn um die Pol- 
achse beschreibt, so hat nach (15) die Achse a, der Kompabkugel in jedem Zeitpunkt die Richtung 
der Tangente an diese Kreisbahn. Die Winkelgeschwindigkeit w, mit der diese Bahn durchlaufen 
wird, ist durch die Umlaufzeit von 84,4 Minuten nach (14b) aus w Tp = 27 mit 


oo Ve = 0.00124 s7 (15) 


festgelegt, falls der Winkel y zwischen a, und der Polrichtung klein ist. Daher wird aus (12) in 
diesem Fall 


3 a oe sin € 
=—)ox i, 0 = —- 
Gg =) X ay Bon sin 
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und daraus weiter ir 
und 


sin € Om. : 
In sin (15a) 
Sin & Mo 


Denn mit der Erddrehgeschwindigkeit @, = 7.29- 107° s und @ nach (15) wird @w/@, fast genau 
gleich 17. 

Diese Gleichung zeigt, daB schon fiir sin y = 1/17, d. h. fir » = 3937 der Winkel ¢, der sich im 
Raumkompa®8 bei der beschriebenen Umkreisung der Polachse einstellt, den vollen Wert € erreicht, 
den der RaumkompaB bei erdfester und schwingungsfreier Aufstellung am Aquator aufweist. 

Wird dem RaumkompaB bei der Aufstellung auf dem Nordpol eine ebene Pendelschwingung 
aufgepragt, so verbleibt @, und damit 1, in der Schwingungsebene, die fiir den auf der Erde stehenden 
Beobachter die Foucaultsche Drehung zeigt, die am Pol mit der negativen Erddrehgeschwindigkeit 
vor sich geht. Aber der Absolutwert |@| bleibt nicht mehr konstant wie bei der eben betrachteten 
Kreisschwingung, sondern schwankt zwischen dem Wert Null in den Umkehrpunkten und einem 
Maximalwert beim Durchgang durch die Mittellage. Dieser Wert |d3| ist beim ebenen Pendel 
identisch mit , wenn y den Ausschlagwinkel des Pendels bezeichnet (Abb. 3). Ist die Amplitude 1 
dieses Ausschlagwinkels sehr klein, so ist bekanntlich das Bewegungsgesetz fiir das Pendel 


Y ==, Sill'® ti, 


wobei w durch w T = 2a mit der Schwingungsdauer T zusammenhangt. Da diese Formel auch 
fir den Raumkompa8 gelten mu, der ja gleichartige Schwingungen ausfiihren kann, so wird 
nach (12) 


A : sin € 
Os) = OLCOS Ot — On — 15b 
| al 0 0 sin é ( ) 
oder 
sin € wo 
—— = — coswt=—17Tr coswt. (15c) 
sin Ey Wo 


Der Winkel ¢ andert sich hiernach periodisch. Er wird Null fiir y = ») (w t = 77/2), d. h. fiir die 
Umkehrpunkte der Pendelschwingung. Die Amplitude von sin ¢ ist im Verhaltnis von 17 sin ¢)5 = 1 
gréBer als die Amplitude v9 des Pendelausschlags. Mit der Umkehr der Pendelbewegung andert ¢ 
zugleich mit a, der Gleichung (12) entsprechend das Vorzeichen. Die Drallachsen 0, und dD, 
(Abb. 1) schlagen nach der Seite negativer ¢ durch. Dieser Vorgang setzt aber eine genau ebene 
Pendelschwingung voraus. Bei der kleinsten Abweichung wird die Bahn zu einer langgestreckten 
Ellipse, a3 erreicht nicht mehr ganz den Wert Null, die a,-Achse mu8 tangential der elliptischen 
Bahn folgen. Daher mu8 an den Umkehrstellen a, und damit der ganze Kugelkérper des Raum- 
kompasses in verhaltnismaBig kurzer Zeit um nahezu 180° herumschwenken, wobei ziemlich groBe 
Lagermomente an den Prazessionsachsen auftreten miissen. Fiir den Grenziibergang von der 
gestreckten Ellipse zur Doppelgeraden kénnen die abgeleiteten Gleichungen keine zuverlassige 
Loésung geben, weil hierbei ¢ den Wert Null erreicht, der bei der Ableitung der Gleichungen aus- 
zuschlieBen war. Der Nebendrallwert der KompaSkugel mit den Kreiselgehadusen ist in diesem Fall 
nicht mehr zu vernachlassigen. 

Bei Aufstellung des Raumkompasses an einem beliebigen Ort der Erde sind natiirlich ebenfalls 
Schwingungen wie bei einem Pendel denkbar. Auch in diesem Fall gilt die Gleichung (12). Aber 
der Vektor a, setzt sich aus zwei Anteilen zusammen, von denen der Hauptanteil jetzt durch die 
Erddrehung bedingt ist. Wenn die Eigenschwingung angendhert eben ist und in die Ost—West- 
Richtung fallt, d. h. in die Richtung des Erddrehungsanteils von Q3, so hat sie keinen EinfluB auf 
die Nordanzeige des Kompasses, wie man leicht erkennen kann. Aber eine Schwingung in Nord- 
Siid-Richtung wiirde auch bei sehr kleinen Amplituden der a,-Achse schon erhebliche Ausschlage 


e ge 2 gebe € I e€ chwingungsanteil von 3 ist 


[3] senr0 =%Wcosoat, 
wahrend der durch die Erddrehung bedingte Anteil nach (10a) 


ls leraar. = Wo COS J 
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\ 


ist. Daw den 17fachen Wert von @ , hat, wird daher in diesem Fall der Ausschlagwinkel « der 
Nordweisungsachse q, in roher Uberschlagsrechnung durch 
YMoMcosMt 17% 


= cos wt 
Wy COS Po COS Po 


tea 


gegeben. Der Ausschlagwinkel « wird also etwa im Verhaltnis von 17: cos Mo gréBer als die Am- 
plitude ») der Pendelschwingung des Raumkompasses. Natiirlich treten dazu noch e-Schwingungen 
gemaB (12) auf. 

Der Einflu8 der Foucaultschen Drehung der Schwingungsebene wurde bei dieser Uberschlags- 
rechnung ganz unterdriickt. Tatsachlich ist er durchaus beachtlich. Denn wenn die Schwingungs- 
ebene zur Zeit t = 0 genau in die Meridianebene fallt, so ist sie nach Vollendung einer vollen Schwin- 
gung, d.h. fiir t = 2 z/w bereits um den Winkel t w, sin y = 2 a(w /) sin p = (2 2/17) sin p oder 
(360°/17) sin p gedreht, d.h. fiir g = 45° um 15°. In (9) ist natiirlich auch die Foucault-Drehung 
enthalten. Aber die strenge Berechnung der Schwingungen aus den Gleichungen (9) und (12) wird 
sehr kompliziert, auch wenn man sie auf sehr kleine Ausschlage der a,-Achse beschrankt. 


7. SchluBbetrachtung. Die vorangehenden Untersuchungen haben gezeigt, da der RaumkompaB 


: gegen Beschleunigungen véllig unempfindlich ist. Seine Vertikalachse (a,) behalt ihre Richtung 


zum Erdmittelpunkt bei beliebigen Bewegungen des Fahrzeugs, das ihn tragt, unverandert bei, 
und seine horizontale Drallachse a, zeigt dabei genau nach Norden, wenn der Raumkompaf auf 
der Erde feststeht. Im andern Fall stellt sich die als Fahrtfehler bekannte Abweichung schwin- 
gungsfrei ein. 

Aber man mu8 beachten, da der RaumkompaB ein Pendel von 84,4 Minuten Schwingungs- 
dauer darstellt, und daf er deshalb vor allen Momenten geschiitzt werden muB, die Schwingungen 
anregen kénnten. Diese Vorsorge gilt insbesondere fiir den Anlaufvorgang. Wie gezeigt wurde, 
reagiert der Raumkompa8 auch bei sehr kleinen Schwingungen seiner Vertikalachse mit beacht- 
lichen Ausschlagen der Nordweisungsachse und der Drallachsen beider Kreisel. 


(Eingegangen am 23. Marz 1959.) 


Anschrift des Verfassers: Professor Dr.-Ing. Walther Bauersfeld, Heidenheim/Brenz, Georg-Beutier-Strafe 9. 
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Zur Membran- und Biegetheorie der Kreiszylinderschale 
fiir ein nichtlineares Elastizitatsgesetz 


Von S. Tameroglu 


1. Einleitung. In dieser Arbeit soll ein nichtlineares Elastizitatsgesetz fiir kleme Verzerrungen 
in der von H. Kauderer! angegebenen Form auf einige Schalenprobleme angewandt werden. 

Zuerst stellen wir das Elastizitatsgesetz in dem fir Umdrehungsschalen iiblichen Koordinaten- ' 
system dar (Ziff. 2) und leiten sodann die Gleichungen her, mit deren Hilfe die Verzerrungen und | 
die Verschiebungskomponenten bei Membranschalen zu berechnen sind, wenn das statisch be- | 
stimmte Problem der Spannungsermittlung gelést ist (Ziff. 3). Nun wenden wir uns dem Problem | 
der biegesteifen Kreiszylinderschale zu und stellen die fiir drehsymmetrische Belastung maBgebende 
nichtlineare Differentialgleichung fiir die radiale Verschiebungskomponente auf (Ziff. 4). Die — 
Lésung dieser Differentialgleichung gelingt nur naherungsweise mit Hilfe eines Verfahrens der | 
Stérungsrechnung. Es werden die allgemeinen Formeln hergeleitet, die man zur Anwendung dieses | 
Verfahrens benétigt (Ziff. 5). SchlieBlich wird fiir den Fall eines biegesteifen kreiszylindrischen — 
Behalters mit eingespanntem MantelfuB gezeigt, wie man bei der Anwendung dieser Formeln vor- | 


zugehen hat (Ziff. 6). 
2. Das nichtlineare Elastizitatsgesetz fiir kleine Verzerrungen. Wir beniitzen die fir eine Um- | 

drehungsschale iiblichen Winkelkoordinaten m und # (Abb. 1) und erganzen diese durch eine zur ~ 
Schalenmittelflache normale z-Koordinate zu einem orthogonalen (¢, #, z)-System, in welchem die 
Spannungs- und die VerformungsgréBen dargestellt werden sollen. Die nichtlinearen Beziehungen 
zwischen den Spannungen und den Verzerrungen Jauten dann allgemein 

Og = 3 K x(é ) & + 2 Gy(ys) (€, — &) » 

O09 ae 3 K x(é) Eo oe 2 G y(y>) (Ep =F Eo) ? (2 1) 

0, = 3 K x(E) & + 2G (ys) (& — €) » 

Too =G v (ps) Yoo ? Tox =6G y (ys) Voz 2 To, > G y(y>) Pox > 


wobei 
pe Ey téa+ ee _ Oo Kk(so) 
0 3 ae PSK 


5 2 
VS = 5 {2 [ey —e0)® + (eo — 4)? + 20] + Yb + Ybe + Ye} 


(2.2) 


gesetzt ist. 

Hierbei bedeuten die Konstanten K = E/3 (1 — 2 y,) den Kompressionsmodul und G den Schub- 
modul. Die GréBen ¢) und y? sind Invarianten der Volumendehnung und des Tensors der volumen- 
treuen Gestaltsanderung. Die Dehnungsfunktion x(¢) und die Scherungsfunktion y(y?) setzen wir 
als Potenzreihen in diesen Invarianten in der Form an: 


(E>) = 1 + x & + x, &2 4 Mz és +--+, (2.3) 
VY) =1 tye Pot ye vet Vers +: 
Man kann die Gleichungen (2.1) umkehren und erhalt so 
if i 
ae (so) 6) + 5G Elta) Open 
1 
Cea F(s9) oy + 3G (4) (3 — 0») » 
1 (2.4) 


1 
Ez = 3 h(S0) 9 + 2G &(%s) (6, — 0) » 


ints 1 
Yos = © 8h) Toa2 Voz = ie 8(t) Tos» RAY tas 


* H. Kauderer, Nichtlineare Mechanik S, 14, Berlin/Gottingen/Heidelberg 1958. 


ee Pee 
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Die in (2.4) auftretenden GréBen s, und t, sind dimensionslose SpannungsgréBen, die man durch 
die Spannungsinvarianten o, und 7, gema® den Beziehungen 
Op + 69 + oO; T) 


oy 
3 Sa SEL Ia yO acai oo 


09 = 


ree (2.5) 


TAG [7 lly — onl + 0 — 04) + (0, — oI] + tho + toe + th | 


ma bestimmen hat. Die Kompressionsfunktion k(sy) und die Schubfunktion g(t?) sollen sich als 


Potenzreihen in den Spannungsgréfen s, und #2 in der Form 
K(so) = 1+ ky: sy + ky s? + hg s8 +-.-, 


(2.6) 
g(t) = 1+ 82: +2, +e: +--- 


- darstellen lassen. 


Fiir zusammengehérige Werte von y? und #2? bzw. von éj und s, gelten zwischen den Funktionen 


—-y(p?) und g(t?) baw. x(e)) und k(s,) die Gleichungen 


1 1 
yw) = g(t) ? HE) = k(s,) . (2.7) 


Die Koeffizienten xi und y; baw. k; und g;, die in den Potenzreihen (2.3) baw. (2.6) auftreten, 
sind zusammen mit den Moduln K und G durch geeignete Grundversuche experimentell zu be- 
stimmen. 


3. Die Spannungs—Dehnungs-Gleichungen fiir ¢, = 0 und ihre Anwendung auf die Membran- 
theorie der Umdrehungsschalen. Setzen wir die Spannung o, senkrecht zur Mittelflache gleich Null, 
so folgt aus (2.5) 


Og +o 


et) K x(é) (3.1) 
und somit aus dem ersten beiden Gleichungen (2.1) 
1 Op + 09 
Oy = ay (6, Sr 65) a 26 vps) ee —— aK x(Eq) 2 
1 op +o 352 
Og = Fou (6, ae 5) a 2 G y(y>) eo ia Ke 5) ( ) 
Too = FVD) Poo - 
Fiihren wir noch die Abkiirzung? 
1 3 K x(é) —2 Gy(y9) 
v (E, 9) oe 3 K x(é ) as G y(y2) 8) 
ein, so wird 
2 
op = 1 — v(ép, v2) G y(ys) le, + ¥(E9> V8) Ea] » ] 
(3.4) 


09 = Sa G y(wo) [eo + r(Eo> V2) Eg] » 


Too =CY(YS) Yoo - 


Die Gleichungen (3.4) kénnen wir nach den Verzerrungen auflésen; dies gibt 


| 1 [oy — (60 ¥8) 09 
“9 T+ le, ¥8)| 26702) |? 
ee 1 [09 — (es 3) op (3.5) 
os A+ ve, ¥3)| 2G yGpf) |" | 
i—— igs 
Poe ~ Gye) 


1 H, Kauderer, a. a. O. S. 144. 
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Schreiben wir mit Riicksicht auf (2.7) die Abkiirzung (3.3) in der Form ; 
1 3 K g(t) cae k(So) (3.6) 


(sorte) = 93 Kgl) + 6 He) 


und gehen hiermit in die Gleichungen (3.5) ein, so folgt 


1 Op ey V(So5 t?) O99 2 
“9 TE WS. #8) | 26 st) 
oe 1 O9 — (Sp, t@) OD (42 Boll 
(oN = 1 a (Sos 12) 2G g(ts) ’ ( ) 


Tod 
Yoo aa G Bite) ov 


Mit Hilfe dieser Gleichungen lassen sich fiir eine Umdrehungsschale die Verzerrungen nach dem 
hier zugrundegelegten nichtlinearen Elastizitatsgesetz aus den Spannungen berechnen. Wir wenden 


nun diese Gleichungen auf die Membrantheorie der Umdrehungsschalen an. 
Nehmen wir — wie in der linearen Theorie — 


an, daB beim Membranspannungszustand allein 
die drei iiber die Plattendicke konstanten Span- 
nungen Oy, Og und Too auftreten, so ist, da wir ja 
auch weiterhin kleine Verzerrungen voraussetzen, 
das Problem der Membranschale nach wie vor 
statisch bestimmt, und die Spannungen in der 
Membranschale bleiben demnach dieselben wie bei 
einem linearen Elastizitatsgesetz. Die Verzer- 
rungen und die sich aus ihnen ergebenden Ver- 
schiebungskomponenten u, v, w sind jedoch jetzt 
auf Grund der nichtlinearen Beziehungen (3.7) 
zu berechnen. 

Als erstes Beispiel behandeln wir die dreh- 
symmetrische Belastung und Verformung einer 
Umdrehungsschale. Bezeichnen wir die Ver- 
schiebungskomponenten in der Richtung des Breitenkreises mit u, in der Richtung des Meri- 
dians mit v und in der zur Schalenmittelflache senkrechten Richtung mit w, so ist hier u gleich 
Null zu setzen, und es wird! 


d 
dp UEP = M1 TE = ¥(Q) » W =T. +g — U cts @. (3.8) 


Hierbei ist r, der Kriimmungshalbmesser des Meridians und r =r, sing der Kriimmungs- 
halbmesser des Breitenkreises. Geht man mit ¢, und ¢, die jetzt nur von p abhangen, aus (3.7) in 
die erste der Gleichungen (3.8) ein, so enthalt diese nur noch gy als unabhangige Veranderliche und 
1aBt sich somit integrieren. Hiermit ist dann auch w aus der zweiten Gleichung (3.8) berechenbar. 

Fiir den Sonderfall der Kugelschale ist r; = r, = a. Der Ausdruck q(y) wird dann unabhangig 
von (Sp, #2). 

Bei der antimetrischen Belastung (z. B. Windbelastung) hat man von den allgemeinen Be- 
zichungen 


ov 
ap = Tee 2 
ou : F 
aa + v cosp + wsing =r, é,sing , (3.9) 


r, Ou en ov : 

ee 7 U COS P + ae = Te Yoo Sin Y 

zwischen den Verzerrungen Eg Eg und Yo und den Verschiebungskomponenten u, v, w auszugehen, 
in denen jetzt p und @ als unabhangige Veranderliche auftreten. Man integriert diese Differential- 
gleichungen wie im linearen Fall mit Hilfe Fourierscher Reihenentwicklungen. 


1 W. Fligge, Statik und Dynamik der Schalen S. 79, 2. Aufl. Berlin/Goéttingen/Heidelberg 1957. 
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Als zweites Beispiel betrachten wir die Kreiszylinderschale. Wir benutzen jetzt statt 3 die 
Koordinate x (in Achsrichtung) und erhalten so anstelle von (3.7) die Spannungs—Dehnungs- 
Beziehungen: 


ee a Oni V(So, t?) Op A 
=—Ttma) | = ae | 
= ras 1 Op — V(So, t?) Ox 
eT wae| ae a(®)). (3.10) 
Tx 
Yep = SP s(t) | 


Auch hier ist das Vorgehen grundsatzlich dasselbe wie bei linearem Elastizitatsgesetz. 


4. Die Differentialgleichung der biegesteifen Kreiszylinderschale. Wir wenden uns nun dem 
Problem der biegesteifen Kreiszylinderschale zu, die durch gleichmaBig verteilte Randkrafte oder 
Randmomente beansprucht ist. Im Hinblick auf eine spatere Anwendung wollen wir vorerst 
einmal annehmen, es liege der Mantel eines Behalters von konstanter Wandstirke 2h und dem 
mittleren Halbmesser @ vor (Abb. 2), der bis zum Rand mit einer Fliissigkeit vom Raumgewicht 0 
gefiillt seit. 

Kénnte sich der Behaltermantel unbehindert ver- 
formen, so entstiinde in ihm ein Membranspannungs- 
zustand mit den Schnittkraften 


n, =oa(l—x), ey S=r0. 


Da die Mantelverformung durch anschlieBende Bau- 
teile behindert ist, so wird der Membranspannungszu- 
stand gestért, und es treten — mit Riicksicht auf die 
Drehsymmetrie — neben den Schnittkraften n, und n 
die Biegemomente m, und m, sowie die Querkrafte q, 
auf. Die Biegemomente m, sind vernachlassigbar klein. 

Fiir das in Abb. 2 schraffierte Element von der Héhe dx 
und der Breite a dp lauten die Gleichgewichtsbedingungen 
fir die Krafte in radialer Richtung und die Momente 
beziiglich der Ringtangente 


n, dy dx + dq,adp —ga(l—x)dpdx =0, (4.1) 
(q, dx —dm,)a-.dp =0. : 
Aus der zweiten Gleichung folgt 


dq, __ d’mx 
Che EP 
und somit liefert die erste Gleichung (4.1) die Beziehung 
d?mx 1 
5 ag (l/— x). (4.2) 
Hierbei sind n, und m, durch die Integrale 
14 “fe 2 Abb. 2 
= ace in, = o,(1 + =)z ds 4.3) 
os f O,: dz Z i ( a ( 


zu berechnen. Mit z ist der Abstand des Schalenpunktes von der Mittelflache bezeichnet. Die 
Spannungen o, und o, folgen aus den Gleichungen (3.4), in denen der Index durch x zu ersetzen 
ist. Die Berechnung der Schnittkraft n, erfordert die Kenntnis der Spannung o,. Wirden wir 
n, = 0 setzen, so wiirde aus der zweiten Gleichung (3.4) 
Cea ae (Eq. 9) Ey 

und somit aus der ersten Gleichung 

Og =o LL + (Em Yo)] 2 © VY) (4.4) 
folgen. 


1 K. Girkmann, Flachentragwerke, S. 448, 4. Aufl, Wien 1956. 
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Um das Schnittmoment m, berechnen zu kénnen, mu zuerst die Spannung o, bekannt sein. 
Sie folgt aus (3.4) zu 


. 2 
* ~~ 1— v(é, Y 


ay © (vs) [ee + Eo Yo) Fl - (4.5) 


Die Schubspannung t,, verschwindet wegen der Drehsymmetrie. 
Mit den Verschiebungskomponenten u, v, w ergeben sich die folgenden allgemeinen Beziehungen 
fiir die VerzerrungsgréBen!: 


Se au . Ow 
Ox ax? ? 
1 dv Zz O?w w 
= eel a 4.6 
ee a op a(a-+ 2z) Doh | are? { A 
tL Ow 2a aeeae Zz Zz Ow | 
Yor ~ G42 dp | a 0x pina aes | 
die sich (unter Beriicksichtigung der Drehsymmetrie) wegen v = 0, la: 0 und OU ae 0 in die 
C oF. op 
estalt 
__ du d?w etd 
Gee ie igs ee ga (4.7) 


bringen lassen. 
Vernachlassigen wir, entsprechend den bekannten Vereinfachungen bei der Biegetheorie der 
Stabe die Ableitung du/dx, so wird noch einfacher® 
d?w w 


ae ee Pe) lol eG (4.8) 


Hiermit wird aus den Gleichungen (4.4) und (4.5) 


Oars Tey [1 + v9 po)] 2 & y(ys) 

4.9) 
is 2 d*w w ( 

Oye 1 — v(E, 2) GC y(ps) |— Seats ape Veo Yo) 

Als erste Naherung nehmen wir jetzt an, daB die Dehnungs- und die Scherungsfunktion die 


Form 
x6) =1, (yo) =1 + yey (4.10) 
haben mégen. Entwickeln wir dann die Ausdriicke 


1 V(Eq, YR) 
1 — v(&, Yo)’ 1—v(Eq, ye) 


2 
V(Eq, Yo) » 
bis zu den in y, linearen Gliedern, so erhalten wir? 


1 
V(Eq, pe) = Ms (L = Vo) (Mia 2%) V2 Yo> 


Lye ad (Iny,) (lee 7 eee 
1 —v(éq, 3) 1 —% L 3 (1 =2y),) : voy , (4.11) 
(Eo, Y5) Yo 1 (lL = 2) 
1 — v(&, 5) el ee 0) i 3 Yo(1 — %) , yovs| * 


Hierin bedeutet 
y —13K—26 
Wo PD Riese E 


1 K.Girkmann a. a. O. S. 469. 
2 W. Fligge, a.a. O. S. 147. 
3 H. Kauderer a. a. O. S. 149. 
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die Querzahl, die sich im linearen Fall ergeben wiirde. Mit diesen Ausdriicken gehen wir in (4.9) 
und sodann in (4.3) ein. Dies gibt 
+h 


I d 
me = 26 | | +m— 3 +m) (20) r¥] 0+ reese, (2) 
=i 
h 
aids 2 (1 + ) (1 — 2») w d?w 2 dw wy 
Ki Somes | 3 (1 — 1%)? 72 Po fon ae mal 1», i 2 dx? 2 a Tall 

=i 

x (1+ yas) (L + fede. (4.13) 


Um nun auch noch die in diesen Gleichungen auftretende Funktion y? zu bestimmen, ist die 
Dehnung ¢, zu berechnen. Hierbei diirfen wir annehmen, da ¢, im allgemeinen klein sein wird! 
gegeniiber ¢, und ¢,, so daB wir wie in der linearen Theorie 


Vp 


hee are (6,71 5) (4.14) 
setzen diirfen. Mit (4.8) und (4.14) ergeben sich dann nach (2.2) fiir ¢) und y? die Gleichungen 
_ 1—2% w d*w 
ores — 9.) (; ee i: 2 GY) 
i 4 2 [dw \2 zw dw , w \2 
Yo = 5 [he (a) 1 Ne Gag cde eld ce: | : (4.16) 
Hierin sind zur Abkiirzung die Konstanten 
_ 2(1—% + %) = 2475 93) 
Us ema (ey rem etme 8S ay) 


eingefihrt. 

Den Ausdruck fiir yj nach (4.16) setzen wir jetzt in die Gleichungen (4.12) fiir n, und (4.13) 
fiir m, ein und fiihren die Integrationen aus. Mit den so gewonnenen Werten fiir m, und n, gehen 
wir in die Gleichgewichtsbedingung (4.2) ein. Dann ergibt sich nach einer elementaren, aber etwas 
langwierigen Zwischenrechnung die folgende Differentialgleichung fiir w(x): 


d? [d2 dw \2 d?w\3 Ew \5 
oy |e (4; — A, w? + Ay wt) +(e) A, w ates As (a) Ap +A, wt — Ayu 
a : d2w\2 d2w\? 
aie az (Bi w? + B,w*) + (ze) (B; w — By w*) — (ae Bs; w* 
+ B,w + B, w — B, w =9 (I—x), (4.17) 


in der zur Abkiirzung die folgenden Konstanten verwendet worden sind: 


_ 81+ %) (lL — 2%) 726 

Ao 9(1—%)? : 
4G h? 

A= 30 =H)’ 


2h 8GY 
ae 3 a lo N2—- 1) yeerny — Ay (H2—1) | ’ 


4A Qh(2h2—a2) , 2h , 2h3 
A; = eal rr ae aes a (211% — mi) » 
if 2 a*h 2h ; 
Ay = Dh oe + [lim — a +h) (4.18) 
2 he 8G. 
Ash [40.— samp 
8 Ag. h? 
y= AS 
2h 2h 2 h3 8G 2% 
4, = (G5 =< 5 at) | 4 qemiee 
32 Ay yh? ni 
A, = 9 (a2 — h?)3 


1 H, Kauderer a. a. O. S. 145. 
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und 
16G(L+%)?%oM2 [2h 2ah_ , 2h® 
fies 9a a Ry ee NS 
po EEE) 2%) 73 8 AP m1 Ne 
age ol (a? — h?)8 ? 
32 G (1 + %)? Yah? m4 
Dea 27 a? , 
326 (1 +») (L—2%) 92 (98 + 2n}) [2ahQh—a) | 2h | 2h 
Ci x Mle (a@&—h?? la ces) aos le 
(4.19) 
P—2 Z 2a°h 2 h® 
Bee 64 Gil oF DNs %) V3 (2 Phe eae a8 —\ ny No » 
4G(1+%) /h , he 
Pe a ( Fas Melredhy 
BR, a 22 EF 0) ah 
7 9 (a? — h?)? ? 
B, — SAE G + %) A = 200) 72 bh (a* + BY) ni 
e 27 (a? — h2)! : 


Auf der linken Seite von (4.17) ist der Ausdruck in der eckigen Klammer gleich m, und der 
ganze Ausdruck dahinter gleich Ty: 
Setzt man in (4.17) y, = 0, so wird 

Atm Aji Ay Ay An AeA Ae ae 
BeBe By = Ba be Be 0 gi Or 


und wir erhalten dann die bekannte Differentialgleichung bei linearem Elastizitatsgesetz. 


5. Die Randbedingungen und die Lésung der Differentialgleichung. Um die nichtlineare Dif- 
ferentialgleichung (4.17) zu lésen, beniitzen wir ein Verfahren der sukzessiven Approximation. 
Hierzu bringen wir die Gleichung zunachst auf die Form 


dtw 
A, Za t+ Bew—e (l—x) = 
d? [d?w 
‘dx? | dx 


d? 
= £5 (8,18 + Bout) —( 


d2w \2 ; dw d?w\5 
(Ay wt — Ay ret) — (“5 ) Aa + (Ge) te + (Gr) to — Av + 4g | 
d?w 
dx? 


2 d2w\3 
J) (Bye — Byw!) +(Ge) Bout — By w+ By w'. (5.1) 


Eine Naherung nullter Ordnung liefert die inhomogene Differentialgleichung 
dw 
1 dx 


A Dal? =o) x 0 (5.2) 


die man durch Nullsetzen der linken Seite von (5.1) erhalt. Thre Lésung ist mit den Integrations- | 


konstanten ¢; sofort geschlossen angebbar. Sie lautet 


Ww 


o - & ] 
sEn eo —32) 


wobei 


ajo Be _ 30 —%) 


gesetzt ist. In dem Ausdruck fiir den dort auftretenden Koeffizienten B, ist, wie der Vergleich 


mit (4.19) zeigt, h3/3 a? gegeniiber h/a vernachlassigt worden, da h <a ist. Die Lésung (5.3) ist 


zugleich die Lésung der linearen Theorie. 


Fiir die weitere Untersuchung sehen wir von der (nur im Hinblick auf eine spatere Anwendung) | 


gemachten Annahme, da der Behalter mit Flissigkeit gefiillt sei, ab und setzen daher von jetzt 


t e4* (c, cosA x + cy sind x) + e** (cs cos Ax + e, sind x) , (5.3) 


ads ah? (5.4) | 


| 


; 


| 


XXVII. Band 1960S, Tameroglu: Zur Membran- und Biegetheorie der Kreiszylinderschale 379 


an @ gleich Null. Die Lésung (5.3) setzt sich dann aus zwei jeweils mit Abklingfaktoren versehenen 
periodischen Ausdriicken zusammen. Der erste (mit dem Faktor e*) nimmt vom oberen und der 
zweite (mit dem Faktor e—’*) vom unteren Rand des Zylinders aus ab. Da die Héhe des Zylinders 
meistens so groB ist, da die beiden von den Randern ausgehenden Stérungen sich nicht mehr 
gegenseitig beeinflussen, so diirfen wir bei der Aufstellung der Bedingungen fiir den oberen Rand 
cz = c, = 0 und fiir den unteren Rand c, = c, = 0 setzen, so da® sich die Konstantenpaare ¢,, ¢, 
bzw. c3, ¢, unabhangig voneinander berechnen lassen. 

7 Wir nehmen nun an, allein der untere Behilterrand sei belastet, und setzen demgemaB 


tr ea On 
Dann wird aus der Gleichung (5.3) 
ww’ — e='* (6. cos A x +c, Sih A x) (5.5) 


Die Integrationskonstanten sind aus den Randbedingungen zu bestimmen. Fiir gleichmaBig ver- 
teilte Randkrafte R lauten diese 


Py 25 Os faroe GQ) ,o9 > 7 (5.6) 
und fir gleichmafig verteilte Randmomente M 
(m) 5 =M,  -@.)25 = 2. (5.7) 
Wirken R und M zusammen, so muB also 
(m,), 5 =M, (%),.9 =—R (5.8) 
werden. 


Wir nehmen zuerst an, es seien allein die Randkrafte R wirksam. Die Gleichung (5.6) befriedigen 
wir dann naherungsweise dadurch, daB wir fiir das Moment m, den Ausdruck beniitzen, der sich 
ergibt, wenn man in die eckige Klammer von (4.17) anstelle von w die Naherungslésung w®) gemaB 
(5.5) einsetzt. So ergeben sich die folgenden beiden nichtlinearen Gleichungen fiir c, und c,: 


— 2) ¢, [A, — A, c? + A, cf] + A, 4/4 c2 c? + 8 A, AS c3 + 32 Ag AM- cS + Apc? — Age? = 0, 
(5.9) 
(cs + c,) (A, — A, c2? + A, cf — 4 A, 2? c3 oc, — 12 A; 4 c? — 80 A, AS c?2) 223 
+ (ce, — €) (4A? Ay 4 cg — 8A? Ag cg 8 + 12 AA 2c? +3 A,e2—S5 Agct)A =—R. (5.10) 


Mit den hieraus zu bestimmenden Werten von c, und c, setzen wir w) nach (5.5) in die rechte 
Seite der Differentialgleichung (5.1) ein und erhalten so die Differentialgleichung fiir die Naherung 
erster Ordnung 


dw 


a dx* 


+ Bw) = os {(c, sin A x — c, cosA x) 2 A? e—4*[A, e—?4* (c, cos x + c, sind x)? 


— A, (cz cosA x + ¢, sin A x)* e— 447] 

— (cz sind x — c, cos A x)? 4 A4 e—?4” A, (cs cos x + , sin A x)® e—34* 

+ (ec, sin A x — c, cos / x)? 818 e—3/* A. + A, 

X (cg sin A x — c, cos A x)® 32 21° e—54* — A, (cs cosA x + c, sind x)® e—34* 
+ A, (c, cosA x +c, sind x)° e— °4* } — (c, sin A x — e, cos A x) 21? e—4* 

x [By (ce, cosA x + ¢, sind x)* e—44* + B, (cg cosA x + c, sind x)? e—?4*] 
— (cg sind x — c, cos A x)? 414 e—?4* [B, (cz cos x + €, sind x) e—4* 

— B, (c, cosA x + c, sind x)® e—34* + By (cg sin A x — e, cos A x) 

x 81S e—34* (c, cosA x + €, sind x)* e—?4% 


— B, (ec, cosA x + ¢, sin A x)8 e—34* + Bg (cg cos x + & sin A x)> e—54* , 
(5.11) 
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Eine Umordnung der rechten Seite und die Ausmultiplikation der Potenzen der trigonometri- 
schen Funktionen ergibt nach einiger Zwischenrechnug 
A, ase Bw = ey oe 34x1(d, — e24* d,) cos 21 x sin A x — (dg — e?4* d,) cos 2A x cosd x 
+ @ 24% d, cos 4/1 x cosAx—e?4* d,cos4A x sind x 
+ (d, —e24* dg) sin 2A x sind x + (dy + e°74* dio) sin 2A x cosAx 
+ ¢24* d.jsin4A x cosdx—e?4* d,,sin4Ax sind x + (d,g + e74* dj,) 
x sind x — (dy, + e°24* dj,) cos A x] } + e 34 (d,, + e 24% d,,) cosA x sind x 
+ (dig + €°24* dyo) cos 2A x cos A x — (dy, + e?4* dy) sin 21x cos 21% 
— (dag + €24* d,,) sin 2A x sind x + (d+ e°?4* d,.) sin 4/1 x.cosh x 
+ (dy, + e 24* dog) sin 4A x sind x + e?4* doy cos 4A x cosAx 
+ ¢ 24% d,, cos4A x sind x + (dg, + e°74* dz) cosA x 


+ (dss + €?4* dyq) sin A x] . (5.12) 
Die Koeffizienten d; sind hierbei 
dy = (22 Ay Cy — Ay C4 — 05 819 Ag) (ef —e) , 
dy = > (2A? Ag Cy — Ag Cy + 32 ¢y Aq 2°) (ef — cl), 
ihe hors b, Ipc e; BA A) =e), 
ped 


= (27? A, ¢, +- Ages + €, A, 324") (c§ — €3), 

= Fear Ag ¢, + Ag €g + 3 A, 414 — c, Ag 32 A) (c§ — 6 2 ct 4 ch), 
d, = QR Ag Cs — Ag C4 — €, Ay 414 — cg Ag 32 A}) (c§ — 6 c2 c? + cf), 
d, = Cs ¢y (20? Ay cg — Az e, — €, 8/8 A;) , 
da Cp 0, (6h 42) (2A7 As ta-— A, 01H Cpe Ag ane) 


d, = ¢,.c, (84° A, cy — A, ¢,— 2A? A, ¢,) , 
dig == Cg Ca (C2 + €2) (2 A2 Ag cy + Ag cg + €, Ay 32 A1) , 
; | 
dy = =F o3 &4 (65 — c?) (2A? Ag cy + Ag cg + €3 Ay 4A* — cy Ay 32 A") , 
1 
dy = > C3 Cy (¢3 — ¢4) (2 A? Ag cg — Ag Cy — Cy Aq 414 — 5 Ag 32.119) 


dyy = (¢8 + 8) (2.48 Ag cg — Ay Cy + 05-4, 829) , 
di =+(3 4+ €2)? (96 cs.4,4!° — 6 Ay c322 + 3 A, co, —c, Ay 4A), 
dy, = (0 + €8) (2A? Ay cy + Ay Cy + & As 82%), 

=> + c2)? (96 c, Ag A}° — 6 Ay cy 2 — 3 Ay cg + ¢, A, 44), 


dy, = 22 cy ¢, (By cy + 22? By ¢,) — > By ¢ (c3 — cf) » 


1 
dig => (cf — ef) cy By + 65 6, (68 + ¢2) (By ey — By 05.2) 22, 
‘ (5.13) | 
dyy = 22 cg Cy (By ¢3 — 2 2? By cy) —> By cg (c3 — ct) , | 


1 
doy = > (¢3 — ¢4) ¢3 Bg + ¢3 ¢4 (¢3 + €%) (By cs + By ey 2A?) 2?, 


d,, = 2? (c2 — c?) (B, ce; — 2 A? Bs c,) + c2c, B 
1 
dy. = ah (c; — ef) (By cs + By cy 2A?) — cf cy (c2 + 2) Bg 


oS 4 
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dy, = A? (c} — c?) (B, Cy 2AAB, €5) 2 eres Be, 


dy, = +P (c3 — ef) (Bye, — By cs 22?) —¢, 2 (ci + cj) Bg, 

dys = 7 c5 ¢, (c2 — c3) (cy By + ¢4 Bs 89), 

dy = + [es ¢4 (ck — cl) (¢y By + cy By 848) —? (By cy + By 04242) (ck —6 c8 ct +8), 
dyy = 7 €3 ¢4 (C3 —c2) (cy By — eg B; 818) , 

dy, = = [5 ¢4 (C} — 2) (¢y By — cy By 8 48) —A? (By cy — By cy 2.22) (cf —6 c2 ct +-e8)], 


“yy = Cs C4 (03 — c2) (By ¢y + By 04222) 22 + + (6g By + Cy Bs 828) (cf — 6 2 c? + e2), 


dzq = Cg C4 (2 — c?) (By cy — By, cg 2 A?) 22 +(e, B, — eg B,; 8/°) (cf — 6 c? c? + cf), 
dey = —> By cy (c3 + 2), 

dg => (ct + €2)? (3 ¢s Bp —8 &4 By 9) , 

deg = —> B, c, (2 + ¢2), 


1 
. 3 + 2)? (3 c, Bg + 8c; B; AS) . 

Ersetzen wir jetzt noch die in (5.12) auftretenden Produkte der trigonometrischen Funktionen 
durch die entsprechenden Summen oder Differenzen, so erhalten wir die Differentialgleichung fiir 
w) in der Form 

) 
A, a8 + B,w) = e34* S” [(D, + e24*G,) sinnAx + (M, + e?4* N,) cosnAx], (5.14) 
n = 1,3,5 
wobei die folgenden Abkiirzungen gebraucht werden: 


1 
D, == [(—8 d, — 6d, + 6d, + 8d, + 16 d,; — 12 d,;) A? — d,, — dy + 2 dys] , 


2 


¢,= ; [(24 d, + 10 d, — 10 dg + 24 dy + 48 dy, — 20 dg) 1? — dyg — doy + 2 dy] , 


o 
1 
D, = = [(— 18 d, — 18 dz) 2? + dy, — dy + dg, 
G3 pele 16 d, + 30 d, + 30d; + 16 d, + 30d, + 16 dy + 16 dy, — 30 d,,) 1? + dig 
— d,, Ee dog — zp 5 
De = <2. 
1 
Gs = 5 [(50 ds + 50 dy)? + dog + doo » (5.15) 
M, = + [(6d,—8 4, + 8d, — 6 d, — 12 dys — 16 djs) 1? + dig — dog + 2 dai], 
N, =F 10 d, + 24d, — 24 d, — 10 d,, — 20 diy — 48 dyg) A? + dog — dog + 2 dys » 
M, = [(— 18 d, — 18 dy) 22 + dyy + dap + Ay 


N, = 5 [(30 dy + 16 d, +16 d, — 30 dg + 16 dg — 30 dy — 30 d, — 16 dy.) ? + dyy 


an doy aia dys ae dog </ 
Mm; 
Wy, = eo [(50 d, — 50 dy) A? — dog + dag] - 
26 
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Die Lésung des homogenen Teils der Differentialgleichung (5.14) ist bekannt. Ein partikulares 
Integral der inhomogenen Gleichung kénnen wir in der Form 


wh) = 944 SS {(D,, + 017** G,,) sin Ae (Mo, e 24% N),) cos nd x] (5.16) 


n= 1,3,5 
ansetzen. Die Koeffizienten’Dyn, Gon, Mon und No, lassen sich bestimmen, wenn man mit (5.16) 
in die Gleichung (5.14) eingeht. Fiir n = 1, 3 oder 5 ergibt sich hier allgemein 
A, A* { Don (81 — 54 n? + 4) + 12 n (9 —n?) My, + &74* (Go, (625 — 150 n? + nt) 

+ 20 n (25 — n?) No, ]} sin nA x 4+ A,A* { M,,, (81 — 54 n? + n4) 

—12nD,, (9— n*) + &?4* [No, (625 — 150 n? + n4) 

— 20 n (25 — n?) Gy,]} cosnAx + By (Don + & ?4*Gy,) sin nd x 

+ Be (Mon + & 24* Non) cos nd x 

= (D, + ¢4*G,) sinndAx + (M, + €74* N,) cosndAx, 


und hieraus folgen die linearen Gleichungssysteme 


Don [A144 (81 — 54 n? + n*) + B,] + [12 4,/4 n (9 — n?)] Mon =D,» 


— D,,, [12 A, A* n (9 — n*)] 4+ [A,A4 (81 — 54 n? + n4) + B,] My, = M,, 
Gon [A, A* (625 — 150 n? + n*) + B,] + [20 A, A* n (25 — n’)] No, = G, > (5.17) 
— Gon [20 A, A* n (25 — n’)] + [A, A4 (625 — 150 n? + n*4) + Be] Non =N,.- 
Aus diesen berechnen sich die unbekannten Koeffizienten zu 
Dn 12 A, A4n (9 —n?) 
Piet M, 4,44 (81—54n2-+ n4) + B, 
On | A, At (81 — 54 nn? + nt) + Beg 12 A, A4n (9 —n?)| ” 
—12 A; dn (9—n®) A, A‘ (81 —54n? +n) + By 
A, i! (81 — 54n? + n‘) -- Bs Dn, 
M —12(9—n’)n A, M, 
on —~ 7 4, (61 —54n' +n) +B, 124,4nQ—n)| ’ 
—12A,A4n(9—n?) A, a4 (81 —54n + nt) + Bg 
(5.18) 
Gn 20 A, An (25 — n?) 
(us Nn A, A* (625 — 150 n? + n*) + Bs 
on = 74, At (625 — 150 n? - n4) + By 20 A, din (25——n®) |’ 
— 20 A, M4n(25—n®) A, At (625 — 150 n? + nt) + B, 
A, M4 (625—150 n? +n!) 4+ Be Ga 
Nee .| = 20 A, At n (25 — n?) Nn 
on A, Mf (625 — 150 n? + n4) | Be 20 A, Ain (25 —n®)| 
— 20 A; An (25 — n2) A, At (625 — 150 n? + nt) Bg 


Haben wir diese Koeffizienten bestimmt, so lautet die allgemeine Lésung der Differentialgleichung 


(5.14) 
wi!) = e4* (VcosAx + Wsind x) + 734% » ([(Don + & 24*G,,) sinnAx 


n=1,3,5 


ot (Men es 4S INg, cos mya (5.19) 


Die Integrationskonstanten V und W sind wiederum aus den Randbedingungen (5.6) zu er- 
mitteln. Setzen wir w"!) aus (5.19) in die eckige Klammer aus (4.17) ein und benutzen den so ge- 
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wonnenen Naherungsausdruck fiir m,, so liefern die Bedingungen (5.6) die folgenden beiden nicht- 
Iinearen Gleichungen fiir V und W: 


{2722 W—7 _&, [ln?—9) Mon + 6 2 Don + (n® — 25) Noy, + 10 2 Gon} 


x {hi—AnIV +S (Mon + Nonll* + 4g (V+ 3 (Mon + Nowll¥ 


n= 1,3,5 


A[V+t > (M+ NP 2AW +2 > -[(n?—9) M,, +6 Dor 
330%. Wi 3.0 


n= 1,3, 


+ (n? — 25) No, + 10nG,,)2 + 4,22 W +42 > [(n?—9) My, +60Don 


n=1,3,5 


+ (n? — 25) Non + 10nGy,]} + A, 2RW +R SY [(n2—9) My, + 6nDon 


n=1,3,5 


=r (n? ie 25) Non te 10n Gonl}* = A,[V 35 ws _ (Mon se Non) |? 


mete bi 4 ee on 1; Noll =O, (5.20) 
{28 (V+W)+2 > [9 Mon (n®?—3) —n Dp, (n?— 27) + nGpo, (75 —n’) 
n= 1,355 = 
pa eT lar as 2 (Mon 1 Nol 
n= 1,9; 
$4,1V+ 3 (Mon + Not— 24,2 PV +R S [(—9) Mon +6nDow 
n=1,3,5 n= 1,3, 


+ (n#—25) Now +10mGoul][V + 3 (Mon + Non) 


=3A4,[22W +2 > [(n?—9) My. + 6n Don + (n?— 25) NM, +10nG, J 7 


Te 85350: 


—5A,[2AW+H Dd [(n?—9) Mon + 67 Don (n? — 25) + Non + 10n Gy ,]]4 
w= 1,3,5 


ee NEV cee eon, «2 3’ Mg oe Gh oa Noel 


1,3,5 


x {2ZRW +R [(n? — 9) My, + 6 n Don + (n? — 25) Non +10 nGon]] 


7 1,350 


X[24,{1V+ 2 (Mon + Non)]—4A43[V + X (Mon + Non)]" 
n= 1,350 n= 1,4, 


+34,[V+ > (Mo.+ No.) P22RVW +2 > (7? —9) Mo, 
n=1,3,5 n 


=1,3,5 


BP6 n 5,, © (n* — 25) No, 4-10.” Cos]? [V7 4 > (Mon + Non) /[3 4,54, 


n=1,3,5 


XV 3S (Mon + Now} =— RB. (5.21) 


Hiermit ist die Naherungslésung erster Ordnung fiir die Differentialgleichung (5.1) unter den 
Randbedingungen (5.6) gefunden. Das Verfahren der sukzessiven Naherung 1aft sich in ent- 
sprechender Weise wiederholen und erméglicht so grundsatzlich die Berechnung der Schnittkrafte 
und der Schnittmomente m,, n,, q, sowie der Spannungsverteilung mit beliebiger Genauigkeit. 

Hat man w(x) hinreichend genau ermittelt, so folgen m,, n, und q, aus den Gleichungen 


d?w\2 d*w\3 d?w\5 
a d’w (4, — A, w? + As wt) + (Ge) A, wt — (5 4s — (Zs) A, 


emmy 
+ A,w?— A,w'*, 
aie 5.22 
UP Ne ee ( ) 
2 Ww 
ny = “8 a (By wu? + By!) + (Fs) a (By w — By w®) — a(S B, w? + B,aw 
Xe 


+ B,aw*?— Baw’. 
26* 
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Ist der Rand allein durch Momente M gleichmaBig belastet, so lassen sich fiir die Rechnung 
dieselben Formeln benutzen. An die Stelle der Randbedingungen (5.7) tritt dann (5.6), und dies 
hat zur Folge, daB auf den rechten Seiten der Gleichungen (5.9) und (5.20) das Moment M und auf 
den rechten Seiten von (5.10) und (5.21) der Wert Null erscheint. 

Entsprechend sind bei gleichzeitiger Wirkung von R und M die rechten Seiten von (5.9) und 


(5.20) gleich M und diejenigen von (5.10) und (5.21) gleich — R. 


6. Der kreiszylindrische Fliissigkeitsbehalter mit eingespanntem Mantelfu8. Das in Ziffer 5 dar- 
gestellte Berechnungsverfahren laBt sich zur Ermittlung der Schnittkraft und des Schnittmomentes 
eines Fliissigkeitsbehalters mit eingespanntem unteren Rand verwenden. Hierbei ist, wie bei 
allen nichtlinearen Problemen, zu beachten, daB das Superpositionsgesetz nicht mehr gilt. 

Im Fall des Hookeschen Gesetzes lieBe sich aus der hier bekannten allgemeinen Lésung fiir 
w(x) auf Grund der Einspannbedingungen am MantelfuB 

dw 
fo(x) pio = 0, wel a (6.1) 


die Schnittkraft R und das Schnittmoment M sofort berechnen. 
Im nichtlinearen Fall kénnen wir nun folgendermaBen vorgehen: Wir setzen die nullte Naherungs- 


lésung 


w\")(x) == ; = o (lL— x) + e%* (cs cosA x + cy sind x) (6.2) 


in die Randbedingungen (6.1) ein und bestimmen hieraus die Konstanten c, und c,. Die so ge- 
wonnene Lésung (6.2) setzen wir in die rechte Seite der Differentialgleichung (5.1) ein und be- 
stimmen deren allgemeine Lisung. Sie enthalt wieder zwei Integrationskonstanten V und W, die 
aus den Randbedingungen (6.1) zu bestimmen sind. Mit der hiermit gewonnenen Naherung erster 
Ordnung fiir w geht man nun in die Gleichungen (5.8) ein, aus denen sich jetzt M und R berechnen 
lassen. Die Zwischenrechnungen sind zwar elementar, aber ziemlich umstandlich. Sie sollen daher 
hier nicht wiedergegeben werden. 


(Eingegangen am 26. Februar 1959.) 


Anschrift des Verfassers: Dog.-Dr. Miih. S. Tameroglu, Istanbul (Tiirkei) Giizelbahce, 
Yeni Kardesler Sok. Refet Apt. 1/7 Nisantas. 
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Bending of Simply Supported Rectangular Plates With 
Clamped Portions Along Arbitrary Sections of the Edges 


By M. Kurata 


1. Introduction. It seems that problems on the bending of rectangular thin plates have been 
worked out for almost all combinations of edge conditions. Nevertheless, problems such as stipulated 
in the title have been left comparatively untouched. In this paper, a method of solving such 
problems is presented in as general a form as possible in an attempt to fill up this gap. The theo- 
retical solution is obtained by means of determining the resisting moments introduced along the 
clamped portions of the edge which shall just cancel the slope of the surface of a plate. 

Numerical results shall be given for three square plates subjected to a uniformly distributed 
load. The first example is a square plate having a clamped portion in the middle of one edge; 
the second example is one having the same portion as before on every edge; and the third example 
is a plate clamped in the neighborhood of every corner. Finally, the experimental results are briefly 
referred to and it is shown that an agreement is found between theory and experimental data. 


2. The simply supported plate bent by a single moment acting upon the edge. The expressions of 
the deflection surface into which a simply supported plate is bent by a single moment applied 
to a given point of the edge are introduced here for latter convenience. Referring to Fig. 1, the 


Y 


0 


Fig. 1. Plate under a unit load and system of coordinates. 


deflection function w of the plate subjected to a unit concentrated load applied to a given point 
P(xo, Yo) is expressed in the following manner: Denoting w by w, for x) < x , we have 


te b2 co . . : 
wy a LAIR Tig SULT {sinh n x (A —&) [(1 + nad coth nz) sinh nx & — 


~ Dx8 — n3 sinh nwa 
—nzé, cosh nx &,] —na (A —&) cosh naz (A —&) sinhnz&\, (2:1) 
and denoting w by w, for x) > x, 
ey b2 co . . ’ 
W. =a pS Ais ee at {sinh na (A —&) [1 + naz coth nz A) sinh na & 
—nzé cosh nz é] —na (A—&,) cosh na (A —&,) sinh nz &} . (2.2) 
In the above expressions, it is noted that 
gee ee, pel: — Yo ae 
Se U lace Ae: aa ee St Saat i=, 


and D is the flexural rigidity of plate. 

Now, using the expression (2.1), the deflection function w, of the plate bent by a single bending 
moment M, what is applied at a given point y = yo on the edge line x = 0 and in such direction 
that the axis of the moment M, coincides with the y axis can be given by the following operation:? 
ain, 


y= Mal. 


x 


%=9 


1 K.Girkmann, Flachentragwerke, p. 195, Wien 1954. 
2 4, Nadai, Elastische Platten, p. 162, Berlin 1925, 
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Table 1 
Case fy(*) fi (x) and their Fouriers representations 
| NI x 
h 
(1) a Sas b x nm 
= ae 
jo) = an haa an Baa ee cosh b (a—x) 
nsin b b 
4 a? m 5 DURES 
mb > 5 as een a 
(1,2, 3, .«.) (m es =) 
1 nza cock : n | 
/ 4 5 
filx) = - ——— + cosh =~ (a— x) — nc sinh = (a—x) 
nsinh sinh 
SS ai Tua nt x x MGS 
Q) file) = 2 footh "54 sin 2F* — = cosh 5 


F n 
(3) _ “i l sinh ; 
Flx) + f(x) = Bac epealeg pe Apt 
n cosh pale Z cosh sthded b \2 
26 26 
_ 8a m mmx 
m7 b oF a2 2 2 we a 
(135) (m ae | 
It x 
fala) + f(x) = : nna OS tha ee 
cosh 26 pe 5 ae (Ny 
26 ery Ts 
nm 
— cosh alg x) 
s ful) —File) = —* ee 
n\*) ~ In\®) = 7) -++ — cosh Ra A 
_, nwa he Se aC a b (3 
h 
nsin 2b sinh om 
_8@ m mmx 
abs = . ae ; sin - 
(G4e yee | 
NnUXx 
n@) Fale) = —— Sa eo 9 
nsinh “7 @ ee inh GS b 77) 
2b 2b 
a UII oe nw /a 
b sinh b E *| 


- 
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The deflection function w, produced by a single moment of the same kind M, applied at a point 
Y¥ = Yo on the edge x = a is similarly given by using the expression (2.2) as follows: 


w, = : 
uv 7a OX Jing 


The same process can also be applied in obtaining the deflection function produced by a single 
moment of the same kind M, or M, which is acted upon at the y = 0 or y = b, respectively. For 
reference, the results through actual operations will now be enumerated below. 


1) The deflection and the slope displacements produced by the moment M,, are 


LAE  WEYo 5 Wty: : 
=D Dy on , sn; Srl) » 


ow M. e oItae 5. WOE , 
oe Re > sin ==" sin Sul) » t (2.3) 
= = Fr > nsin 22% cos nee. f,(2) . | 

n 


2) The deflection and the slope displacements produced by the moment My, are 


_ My 5s WORE) ona ee } 
0, => > sin sin = bal) 
n 


a 
dwy aS My - NUX% . NUXX / 2.4 
Diese care een tue, an 
dwy My - NIX OEE OE 
Dee Fy gos - Bul) + | 


The function f,(x), f(x) contained in the expressions (2.3) are those shown in the item of Case (1) 
in Table 1 together with its Fourier series expression. 

Next, the deflection and the slope displacement produced by the moment M; are also given 
in the same forms as expression (2.3). In this case, however, it should be remarked that M,, in (2.3) 
is changed to M; and the functions f,(x), f,(x) shown in the item of Case (2) in Table 1, are 
used instead of f,,(x), f(x), respectively. 

In the case of symmetrical bending when the plate is acted upon simultaneously by both mo- 
ments M, and M, whose magnitudes are equal to each other, the deflection and the slope dis- 
placement are given by expressions which are replaced respectively by f,(x) + fr(x) and fn(x) 
+ fi.(x) shown in the item of Case (3) in Table 1, instead of f,(x) and f,(x) in (2.3). 

In the case of antisymmetrical bending of the plate subjected to moments M,, und — M, of the 
same magnitude, the expressions in (2.3) become applicable by replacing f,(x) and f,(x) by f,(x) 
— f,(x) and f.(x) — f,(x) shown in the item of Case (4) in the table. 

The expressions by moment M, or M, are shown in (2.4) ,which are obtained from the ex- 
pression in (2.3) by replacing M, or M; by M, or M, and the transpositions (a, b) and (x, y). 
Therefore, the functions g,,(y), g,(y) in (2.4) are, corresponding to each of the cases, those produced 
from f,,(x), f(x) and so on tabulated in Table 1 by similar rewriting. The Fourier series expressions 
of f,(x) shown in Table 1 can also be written briefly as follows: 


2 
fila) = 35 een = Se Sa) 2:5) 
te BPA (m + — n*) Z 


in which the coefficient k,, signifies values tabulated in Table 2 corresponding to each case. 


Table 2 
Case | 1 | 2 | 3 | a 
ht Pte eee f1-cy 
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3. The slope displacement at the edge produced by the single moments acting simultaneously 
on several points on the edge line. The following nomenclature is used with respect to a single 
moment acting on a given point of the edge as below: For instance, M.,; means a single moment 
which is acting upon point y = v; on the edge line x = 0 and in a plane perpendicular to the y-axis. 


M,, means a same moment on the opposite edge x = a. Similarly, M,, means the moment on the | 


A U 
edge y =0 and M,, means that on the opposite edge y = b (see Fig. 2) . In the next place, let us, 


for instance, denote the slope displacement at the edge x = 0 caused by the above-mentioned 


moments as follows: @,.; is the slope displacement caused by M,,, On 7, is that by M;,, Oxy, DYM 


M,, and On9; by M,,- 


Fig. 2. Notations concerning edge moments, 


Using these notations, the resultant angle displacement 0, _ produced by the combined action 
of these moments can be written in the following form: 


a B y 6 
O,=0 = 3 Oxn, + 3 On3, + dy) Gey +d Ong, (3.1) 


in which each summation sign on the right hand side means the sum corresponding to the number 
of single moments on each edge. Similarly, the following expressions are obtained with respect 


to the resultant angle displacements: — 0,_, at the edge x =a and those at the edge y = 0 
and y = b, 


a B y 6 

Oz=0 = ee, Oz., am > O35, Te BNE) a 2 Oz, ’ (3.2) 
a B y 6 

Oy 0 = F Oye, + Oya, + Y yy, +S Oyg,> (3.3) 
a B Y 6 

Oy 20 = 3/954, + 2 Ops, + DOs, +d O5,, - (3.4) 


Each term on the right hand side of the above expressions may be written down by using (2.3) 


or (2.4). Taking each term of (3.1), for example, we have by rewriting the subscript of w,, M, 
by x; and taking as y, = v; 


rat i re eer a oe Ue UCT ane 
oF | Pee ie > sin —— sin — =f Os (3.5) 

Ow; Mz. nv. eee a 

“ 2 = u c : Dairy fk 
0, 3; aX) |e=9 4 Da SS eth ate b Jn(0) (3:6) 


in which f,(0) is that obtained by taking x = 0 in the expression f.(x) of Case (2) in Table 1. Simi 
larly, rewriting the subscript y used in (2.4) by y; and taking as x) = u; or rewriting that by y 
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and taking as x, = u,, the following respective expressions are obtained: 
Owy. My, nu 
Sauk Ni oes Yi . t 
eu, =| Se leno = Da Oy tM Bal)» cet 
| Owy; Mj; _ neu, _ 
0) =| Be lena = Da Dy POG BaD) CB: 


where g,(y) is that obtained from f,(x) of Case (2), Table 1, by such transpositions as already stated. 
Hereafter, it will generally be considered that the top bar is attached above the notations repre- 
senting the quantities relating to the opposite sides against the coordinate axes. 

Next, for each term of (3.2), similar derivation from (2.3) or (2.4) also gives 


15 = |e bee De Date tin FE Sle) 9) 
Ie fo: ae = De tin or tain ne f(a) , . (3.10) 
Gy) ae an), (3.11) 
Ce) Oe ites a) (3.12) 


n 


Since the meaning of (3.3) and (3.4) is quite analogous respectively to those of (3.1) and (3.2) with 
the exception of the difference of the concerned sides of the plate, each term in the right hand 
side of such expressions may also be derived from (3.5) to (3.8) and (3.9) to (3.12) by making the 


following transpositions: 
(x), Une), (May My) (3.13) 
(a, b) ? (u;, v;) ? (M;., M,,) . 
Accordingly, the rewriting of the functions must be carried out by such transpositions as 
frl0)—> gn(0), — fn(0) > Bn(0),—Ealy) > Fale)» Enly) > fal) - 


_ 4, The slope displacement at the edge produced by the moments distributed in several sections 
along the edge line. Next, let us regard the single moment M M 2, M,,, M,, respectively, as 


the intensities of the distributed moments in the intervals h,, h;, |, 1;. Then the expressions for the 


total slope displacement at the edge produced by such distributed moments can be obtained by 
integrating the right hand side of Eqs. (3.5) to (3.8), Eqs. (3.9) to (3.12) and so on through their 
own intervals and substituting them into expressions (3.1) to (3.4). The actual operation gives 


a B ea Teer 
CO, =| DXi sinnzy + DD, XnihelO) sin nn] 


7 


ay as a Yui” 8nl) + 2 D>) Yn 70) (4.1) 


z 


de Dy Pa) Ga Deh Sn) =F 25 Pa Yo 1)e nat (4.2) 
Oyo OS Y,,, (0) sin nvé + YD ki sinns| 

Tee Ds Dy Xai mfal€) + Ss dX nh| (4.3) 
Oyat =7|> pai Ynign(0) sin nc é + > 2s Yn Bn(b) sin nz | 

+ nS oS Xni (— 0 fal6) + a Dy i (—1)" nit| (4.4) 
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where 
x y 
E = aa p) 1) at b ? 
b/ b; & 
nav. _— Z na v; =a (4 5) 
. i =) af a 2 
A = | M,, sin . dv; 3 1G, — | Me sin b dv; 5 { 
b; b; 
yy —_ 
a, ay as 
nN 7 u. = eh Ue eee 
. v — . 
NG | M,, sin — du; . Yo | M,, sin — du; 
a; Gj 


and a, B, y, 0 mean the number of sections to which the moments are distributed. The upper and 
the lower limits of the above integrations represent the end points of the moment distribution, 
that is é f 5 

36h: a — Oh b; —b, = 1,, ,—b; =1,. 


Now, assume that such distributed moments may be taken in the form of the following trigono- 
metrical series: 


seu. 4s m 

Me 2 E, ; sin Ze 1 Leg oe sin Te ia 
SIX. SMX 
| Po > F,; sin haa M,, = 2 Be sin 


in which : : 
0<y,<4,, 0<y,<1,, Ory eon. 0<— 4, <Ohe. 
Substituting the series (4.6) into expression (4.5), we have 


f = rm 2 f 
Xn =—t >) 5 E,;D,, (6,1), Xi =— 2 D's EG, (6,1), 


: : (4.7) 
l; = i; = ~ 
Ne a >) 8 Fei Pun (a, hi) Se ae >) s Fy: Dn (a, hy) 


where the symbols @,,, used in the above, e. g. ®,,,(b, J;), represents the function of the following 
type: 


(— 1) sin ee sin 


P, n(b, Ii) = 


(4.8) 


From this, it will be considered that @,,,(b,1;) represents the function obtained from expression 


(4.8) by replacing b;, b,, I,, by 0j, 6; 1;. Similarly, ®,,(a, h,) is obtained from expression (4.8), 
exchanging b;, b;, 1; and b, respectively, for a., a;,h;, and a. Moreover, ®,,,(a, h;) means the function 
obtained by taking aj, a;, h; instead of a;, a;, h; in the above function. Hereafter, the symbol ® 
will be used as the expression for the function of type (4.8), taking into account that the use of the | 
other notations for contained variables or suffixes means the replacement of them in the corre- | 
sponding places of expression (4.8). 


5. Modification of the representable interval by a trigonometrical series. Each of Eqs. (4.1) to (4.4) | 
introduced in the previous section represents slope displacement along the full length of the edge | 
on one side. Expression (4.1), for example, represents the slope displacement at any point on the | 
edge x = 0. Now, let us consider for the sake of laiter convenience, the method by which those | 


expressions become representable only in the interval of the distribution of each moment. Take 
as an example the following series: 


H(x,) = > A,, sin “nat 
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which can represent some given function in the interval 0 < x, < a. Expanding again this series 
in Fourier series which is representable in the subinterval ay < x < aj (or 0< x, < h), we obtain 


A. sin “24 = aS m sin “7 Si ®,,,(a, h) , (5.1) 


in which 
(— 1)" si Yes eine 
@ n(@, h) = = = i : (5.2) 
Pe ye 
re 


The left hand side of Eq. (5.1) is an expression of function H(x,) in the interval 0 < x, < a and the 
right hand side stands for the same function only in the interval a,< x,< a, or O< ay<h. 
Making use of relation (5.1), expression (4.1) can therefore be rewritten into the expression which 
holds ealy 4 ind< x< I; (bj —b; =1,). First considering the above, the first term on the right 
hand side of expression (4. 1) is transformed into 


Dd; Xniful0) )sin "EX = — TD mein” mT 13 El) Pryn(b, Uy) « 


The same process is applicable also to the second term. Concerning the third term, if relation (5.1) 
is applied to the series form of g,,(7) (see pee we obtain 


» Yume, (2) = 22 > ys Pee: zi sin ™2Y 


m+ 


ae Dy Tye eh ) (k, = 1). 
if (t+ 3 


The forth term also can be rewritten in the same way as above, provided that k,, = (— 1)™—! is 
taken to make the series representation of g,(y/b). Eventually, we find that expression (4.1) is 


transformed into 
-B Aaa 
*1Y SPS) Xafl0) Baal) +S) S Xa) Palla 


Oo o=— D pal min 


6 
4b? aa (— I)cet r res 
ee i pe ea Ot) (5.3) 


Similarly, expressions (4.2) to (4.4) are transformed into 


es pay renee YD alo Dy +> > Kee f(a) Din lb, I) 


+n > > aden "n>! (on eo Pn (b, L,) 
SEs SY, Hip a > ee) @,,,(b, i) (0<3,<i) (5.4) 
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6 
O,£9=— tee ai m sin a ae Sy) Sy Yn i8n(0) ®,n(a, Aj) + Ds, 2 Yn i Bn (0) Pan(@, hy) 
peas cpap ten 


a 5 ea Aen } = Bees Lees _®,,,(a, h,) (0 ei Xj a h;) f (5.5) 
= (P44 


6 
O25 = a oy m sin — DS, Dy Ynign(d) ) Dy, ,(a, h;) + > a Yn; gA(b) Py n(a, h;) 


b2 
B — eis 1 r—1 — ee as 
te Hp a Xqi(—1)"n 2) (amen h;) (0<x,<h,). (5.6) 
b2 


Thus, we have found the expressions representing the slope displacement @ only in the section 
of the distribution of a moment. 


6. The conditional equations for partial clamping. The deflection function of a simply supported 
plate under a given load is expressed by a single sine series of x or y as widely known, and then 
the angle displacement g along the edge can also be expressed by a similar series. Now let us 
assume that those angle displacements are given as follows: 


1 SeeTUIELY; 
Px =0 = pao B sin 


il TUS ees 1 Fi oa. nae 
Qy=0=z- > 4, sin, one ¥y == 03 Py=6 =Hz Dy Ansin 


Referring to the procedure explained in the previous section, those slope displacements are trans- | 
formed into the following series which are applicable only in sections of moment distribution, 
such as 0< Vp ;: 


iJ ae 
am A= Os P2=0=5z >, Basin, On be 0) 3 | 
H (6.1) | 


on =P 


Tee SY msin a 3, Be Ondo 1) (sey (6.2) 
Px =a = — pesmi” wee Pr rlb, O< ¥j<h), (6.3). 
(a0 =— pos > m sin — yoy ae Pm n(a, (0<2,<h)), (6.4) 
Py=b = 2S min so ES) An Palas ty) — (0< 3) < hy). (6.5) 


If resisting moments on the side edges can be eons to cancel the slope displacement @, _ in the 
interval b;< y< bj (Le. 0<y,< 1 ;) caused by a load, we must have 


pee ena for 0 yy) 


in which 0,9, Y;=0 are given by (5.3) and (6. 2), respectively. Regarding slope displacements | 
on other side edges, the same consideration again holds. That is 


Op a0) Og OO OV i<ohs 
Pr-et+Qraa=0 for 0<¥j<I,, 
yao Opa == aston mene aj<hy, 
Gy=-.+O,-,=0 for O<a,< h, . 


(6.6) | 
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To make each of the equations in (6. 6) into identity, we must equate their respective coefficients to 
zero after the substitution of expressions (5. 3) to (5.6) and (6.2) to (6.5). In this way the following 


equations are obtained, using the expressions [(5.3) to (5.6)] for © which have been substituted from 
expressions . 7): 


for b;< y < bj, (b; — b; =1,), 

DBs Prn(bsl,) - 4 Pa >} fal0) @ Dru nlb, h) 
et Aa ‘She n(b, 1) Dy n(b, 1)) ) 
eA Sia ZI i? 


+ys as wy 2) ae a ae enla, h;) D,, (b,j) | 5 (6.7) 
for bj; << y < bi, (b; Sh =F i 
>» b, Dy n(b, = 2e[ae 1 >) fale) P,,(b, 1;) Dy n(b, f) | 

£ ase DB Dj fala) ®, (bs 1) Pu n(b, i) 


+e dF BE ay ue hi) ®,,(b,1,) 


sys SF Pane Selatan lh an 


for a;< i qj, (a; — a; =h,), 


DA Pandy) = SE [Ss Fas 3) 0) Balto) Pals 
n t a s n 


P,,(a, h;) ®,,,(b, I;) 


ss 
LM 


US sare eee ape aah) Ike (6.9) 
I 2 See ny J 
for a;<*%< Gj, (a, —a, = h;) > 
oS Ae D,,n(a; h,) = . a py s | oan > &0lt) P,,,(a, h;) Dr n(a, hy) 
[Ds Fs 3) EO) Pal fe) Panta By) 
4a? oe eo, 1) Pac, 1) 
45RD Da Bates ny 


=, 


r 


B = = 
Aa? I see nr f 6.10 
+5 = Ys 2) > a ©. (bs 1) D* (a, h;) |. (6.10) 


: : (n+ gem t) 
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The above obtained four sets of equations express the conditions of clamping as explained before. 
The number of equations in each set corresponds to that of the clamped sections, 1. e., Eq. (6.7) 
represents a kinds of equations according to the various values of subscript j. Similarly Eqs. (6.8), 
(6.9) and (6.10) represents respectively / kinds, y kinds and 6 kinds of equations. Moreover the 
system of equations in each kind consists of an infinite number of linear equations corresponding 
to the subscripts m and s. 


If the unknown constants E,,, E,,, F,;, F,; are obtained numerically by solving such a simul- 
taneous equation, each of the moment distributions shown in (4.6) may be determined as a statically 
indeterminate resisting moment. Here it must be noted that the double series contained in those 
equations can generally be transformed into a single series for convenience of computation. Such 
an actual operation in detail will be illustrated in later examples. 


7. Expression of deflection. If resisting moments (4.6) have been determined, substitute ex- 
pressions (4.6) into (2.3) or (2.4) and integrate them respectively through the sections of moment 
distribution, then the deflection of a plate caused by such distributed moments is expressed by 
their sum. For example, let us consider the deflection produced by moment M. a which is given 
by the first expression in (4.6). Integrating the first expression of (2.3), which has been substituted 
from (4.6), gives 


Uy 


P<) 5 WOES? We 10GB (Os FEB) 3 Saree, 
Wa = Dadi dy Baihle) sin “| sin wrigeeoke L dy; 


0 
= — pia D) fale) sin "3% DS? s E,,B,n(b, 1) (7.1) 


where, since b; << yy) < bj, it was taken that y) = b; + y;, 0< y;< 1;. The deflections produced 
by M 2 M,, and M,, are also obtained in the same way. The results are 
I 


; vA ° nN a = 
Ws, = — pia D, Ialt) sin [> D's E,; B,,(b, 1), (7.2) 


h; : 
Wy, rar Dnt >) aly) sun “<= s 1 ®.,(a, h;) 2 (7.3) 

hee te es s | 
So Duk Dy, Bn(y)sin Ds F,; ®, (a; h,). (7.4) 


Denote by w’ the deflection of a simply supported plate which is produced by a given load, 
then the required deflection in the problem concerned is given by the superposition of the above 
obtained deflections in the following form: 


; a B y 6 
ay peach i Poem ae pects Doh (7.5) 


Further, stress resultants such as bending moment and so on can be determined, substituting | 
the above obtained expression (7.5) into well-known formulas. 


8. Illustrative Example 1: Symmetrical bending of a simply supported plate clamped along the | 
middle portion of one edge. Let us suppose that the clamped portion h, is on the edge y = 0 
(see Fig. 3) and the clamping moment distributed on h, can be expressed as 


fe . SUX 
My, = >) F..1 sin ace (8.1) | 


Taking as j =y=1, a =f =6=0 in the conditional e i i 
‘ > quation (6.9) about edge y = ide- 
ration of the fact that there is no clamped portion on other ee a get Meee 


hy, , 
‘ Ds $ Fi >" gi(0) B,n(a, hy) Dy nla, hy) = >) An Pinna, by) . (8.2) 
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If the clamped portion is symmetrically situated on the edge line, i. e., if aj = a —a,, it can be 
expressed that 


; 0, when m-+n_ is an odd number, 
Roeee Ivar G. at 
(— 1)™ sin ——_ — sin ——1 — 
a 


a D . 
—2sin—— ,_ when m-4 n is an even number. 


Fig. 3. Plate clamped partially on one edge. 


Considering the above relationships, expression (4.8) is reduced to 


Te ay 
sin a 
@,,,(a, h,) = — 2 to when m-+n=even, 
nies 
a2 
E (8.3) 
sin : 
@,,,(a, h,) = — 2 a when s-+n=even. 
se St 2 
2 
And on reference to the item of Case (1) in Table 1, we have 
em 1 nazb 1 nb) _ oO; 
g,(0) Shae nab = = ae + cosh am as . (8.4) 
a 
in which 
0, = 20, — coth 2 a, (2.4, coth 2 «, — 1), = ee 


Moreover, assuming a load distributed symmetrically in x-direction, we find that the distri- 
bution of the slope of the deflection surface along the edge y = 0 becomes symmetrical in the same 
direction and then can be expressed by an odd series, 1. e. 

1 a - nwx 
rao = De > A,, sin aii (8.5) 
(1. 3.5...) 
Accordingly, n = odd number. From this we sce that both m and s are odd in reference to (8.3). 
Now for simplicity, let us consider the case of a uniformly distributed load (load intensity: q). The 
deflection of a simply supported plate is given by the following well-known expression, i. e. 


4 qat th Taek a, tanh o, + 2 NI b 
w = —! Se sin te 2 = cosh ™*(y— 3] 


De n? a 2 cosh a, 2 
(153.3...) 
1 nw WN An, He b 
een a (y y sinh . [y 7)| ; (8.6) 
Using this, we find 
sin 2 
__ | dw’ 2qa my a eed ke me aa 
Py =0 ie Lis Dat 2 nt cosh % (siuh On cosh = ‘ o) 


(53.5.2) 
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Then the comparison of both expressions (8.5) and (8.7) gives 
s AGa ii : Cie 2 qa? 

(n= ‘, 35.) wnt cosh a, (sinh pes | Sy ce rent eee Ge 

where 
On = % — tanh «,, (a, tanh a, + 1). 

In the sequel, the conditional equation for the case in which the plate is symmetrically clamped 
along the middle portion of the edge y = 0 and subjected to a uniformly distributed load is deduced 
from Eq. (8.2) into the following form: 


, sin? = ue Ae sine 
On On a 
De a jee 7 Co 3 pe ea (te) =="1,-35-9;-. >.) (Sees 
Ss n m2— =! n2 Fo seca n? 1 n m2 ——_1n2 
(1.3.5...) (1.3.5...) AE ae (1.3.5...) a 


The above equations consist of an infinite set of linear equations whose unknown constants F, , 
are obtainable by solving them simultaneously. 

In the following, such a case of a square plate — the clamped portion is one-third of the length 
of the edge y = 0 — shall be handled as a numerical example. That is 


2 hp et a, i! fi ¢ 
a=); Pitas ee BS eS ay a—2a,=h,. 
And, since 
1 
a, =~ (a—h), 
the relation 
5 LeGEG! 5 UNE n-1 nah 
sin SE oe (Oa) econ ee (aE SSS ees) 


holds. Then, putting h,/a =x, we have for the case when x = m/n 


. naa n—1 nx 
sin i 2 cos —— * pil ee eee 
@ ee (D) 2 Se eS 2 
ling alin = (1) Fate 
m 2 hi 2 m n® m\* 2 n? % 
pips LT ae 1H *—> — Soa, — x? ae 
= a2 a n aces 
5 NU TE 
n—l1 Sin —— m+n 
(22) ) Sie Se ee ee 
( ) ie) ( ) 4m 
4n|— 
n 


Here, considering h,/a=m/n , it is concluded that 


m+n=m{1 +;,)=m +3) =4m 
4 
and then 
m+n 
(1) 32 (ek: 


Thus, the following relations are obtained: 


m i 
when Warm iieas ti 
n a 3 
BEE 
sin i 
a It 
2 poeeal ait ? 
aaa am 
ee 
s h 1 (8.9) 
when — =? = — 
n a O12 
_ nwa, 
in) 
a mM 
2 1 2 As 
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Taking the first twelve coefficients of 5 F s,1 00 the left hand side of Eq. (8.2)’ and calculating each 


of them with due consideration of the above relationships, we obtain the equations shown in 


Table 3. 


Solving of this system of equations yields the following values for the coefficients: 


Ey == 0483 174 Ke, 13 Fis, = — 1.52374 K , 
3 F,, = — 0.611130 K, 15 F,,, =— 1.70191 K, 
5 F5, = — 0.798825 K , 17 F,,, = — 1.88462 K, 
7 Fy, = — 0.982804 K , 19 F.,, = — 2.06582 K, 
9F,, =— 1.16385 K, 21 Peg = 2.23770 K 
LFS — 134416 K, 23 Foy ) = —= 2.42428 K 
where 
3 2 
Kasai 


Referring to the expression (7.5), the deflection function is expressed as follows: 
w= Ww! + wy, (8.10) 


where w’ is given by (8.6) and wy, by (7.3), i.e. 
‘ > s Fea P,,(a, hy) 


(1-3-5...) (1.3.5...) 
5 na x 
=~ pant Dy Far Dd, Bly) Pen(a, hy) sin 22 (8.1) 
(P3050) G:3:5) 


In the above, g,(y) and @,,,(a, h,) are specified as 


sinh = / 
Was — sara + % coh n (1 — 2) : 


8n(Y) ie nsinhnz 


sinh nz b 
( sin = a 
zy 1 
jae ee ss when — zs, 
®,,(a,h4) =} — §—F (8.11)’ 
% S 1 
e as > when a = 730 5 


Calculating the expression (8.11) with the previous values of the coefficients, we can determine 
the deflection by (8.10). Ratios of the deflections of some points to that of the center point are 
described in Table 4 and figured in Fig. 4 by means of contour lines so that the effect of partial 
clamping can be recognized. 

In the foregoing, we used practically the twelve equations for the determination of F,,. But 
an infinite number of equations should be used because F’, , originally meant a term of a sequence. 
Therefore, we must investigate the convergency of such a sequence. F irst, neglect all the other 
terms except the first term on the left hand side of the first equation given in Table 3, then we 
can determine the value of F,, and denote it as the first approximation. Second, take the first 
and the second equations from the table and neglect all the other except the first and the second 

terms on the left hand side as before, then we can obtain the values of F,, and 3 P34, and denote 
them as the second approximation. Afterward, proceeding in a similar way to determine the third, 
the fourth, .. ., etc. approximation, we obtain the results shown in Table 5 (namely, the previously 
obtained values are understood to be the twelveth approximation). Let us calculate the value 
of the deflection at the center point of the plate by using these approximations of coefficients. 
Denote by w, the values of deflection calculated with the use of the first approximation of see 
cients, and by w, the value due to the second approximation and so on, then we obtain the results 
27 
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in Table 6 and their graphical representation Fig. 5. From this, we can see that two significant 
figures at least can be readily decided for the values of deflection. Therefore it may be believed 
that the configuration shown in Fig. 4 is accurate enough from a practical point of view. 


8 
Ss 
= SS 
Sis 
SYS is 
Cunae 
8 
aS} 
ARS 
Sy 
) 8 
Flat fal 
SSESTELTEE 
0 = e © Sssssssgg2l 
i a y-icja ga! 
0.3Wy = We =0.002.96 — 
04M : y=(idJa Y 
05M, = 
0.6NG y=(Me)a 
0.7Ny === Simply supported edge 
is; DIVA Z ——Clamped edge 
Ane Y=Cle)a i Y=(Me)a. 


10M 


Fig. 4. Contour lines of deflection surface for the case when the middle portion of an edge is clamped. 
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Fig. 5. Approximate values of center deflection. (cocte. = De) 5 
I 


Moreover, let us investigate to what extent the condition of partial clamping is satisfied. Using 
(8.10), we see that the slope along the edge y = 0 is given by 
Ow’ 


Ow is 
y =0 oy 


bed Owy 1 
oy 


viene | ey 
in which the first term on the right hand side is calculated by (8.7) and, from (8.11), the second 


term is given by 
= g Cn eeteon Y 
jotoee Dae > sF,, > = ©, (4, h,) sin - 


s 


3 
Sp 


Yu 
oy 


Ee 
1 


and here 
On = na — cothnaz (nz coth nz—l1). 


The actual numerical results, which the use of the twelveth approximation of coefficients gives, 
are presented in Table 7 and Fig. 6. Observing these, we find that the condition of the edge is 
Ih 85 2 


c Or . ° ° x 
fairly well satisfied. Notwithstanding, near the points pias Cis ye where the edge con- 
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Table 7. Slope Distribution along Edge y = 0 [cont = a ) 
I 


D 4 
x/a 
negie 1/12 2/12 3/12 | 3.5/12 | 3.75/12 | 4/12 | 425/22 | 45/12 | 5/12 | 6/12 
2 0.184 0.347 ~ 0.480 | 0.533: 0.556 | 0.577 0.596| 0.612| 0.637] 0.657 
ie) —0.084 | —0.177 | —0.294 | —0.376 | —0.433| — | —0.599 | —0.616 | —0.633 | —0.655 
Sum 0.100} 0.170} 0.186| 0.157| 0.123 || — |—0.003| —0.004| 0.004! 0.002 
e 8 
N “pis 
y 
A 
oe 
= if | 
So) 
2 | 
| 
Le § gsi jasi 8 ae 
S sy eae Eis by al Pz 
2 Tle ae He 35fe | He (se Sle Sie 


Se 4. 25/2 


3 
Fig. 6. Slope distribution along Edge y = 0. (coett = = = ) . 
Ea 


ditions abruptly changes, the above expression become incomputable because the series contained 
in the previous expression at such points do not converge. Now, the clamping moment on the edge 


is given by (4.6). Taking the center points = = = of the clamped portion, for example, we have 
Me = >; Fy. sin a Ee Heyete Pgge = baa 


(4:3:5...) 


This series is converging too slowly to be evaluated, but if Euler’s transformation! of series to im- 
prove the convergency of an alternating series is applied, the value can easily be evaluated. Such 
a transformation is formulated as 


where A"a, means a foreward difference of n-th order. Now, applying this formula to the second 
and the subsequent terms of the previous series, we obtain the following result: 


M,, = — (0,433174 — 0,203 709 + 0,159764 — 0,140400 + 0129316 — 0,122 196 
4 0,117210 — 0,113460 + 0,110860 — 0,108727 + 0,106557 — 0,105403 + ---) K 
Se Ass 17K = 0,203.709 + 4, 00043 945 +1 5, 0,024 581 ++ “a 0,016301 ---- } K 


= — 0,433174 K + (0,101 854 + 0,010986 + 0,003073 + 0,001018 + 0,000374 


+ 0,000148 + 0,000062 +.---) K—- —0,3156 K = — 0,09593 q a?. 
3 1 : 
Similarly, for the point = =F such a series becomes 
My,= >) F,1sin*7 =0,707106 (Fya + Fea — Fea — Fra + Poa + Fua—***) 
(1.3.5...) 


The above series can be regarded as an alternating series by summing every two terms. Then, 
applying the previous formula to this, we find 


M,, == — 0.100q a*. 
1 K. Knopp, Theory and Application of Infinite Series, p. 244, London and Glasgow 1954, 
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At the change points of the edge condition, the above series discontinues to converge ae then 
evaluation becomes impossible. This fact is due to the fact that such points as considere 1 
are singular points where stress concentration occurs. Next, to derive the expression of the poe ing 
moment caused only by the clamping moment, the necessary expressions, from (8.11), are written as 


sar i e ‘ h;) sin 
Gee = Datdy 8 Pat Dm Bul) Byala hy) sin , 
sinh n 7 (2 — ) 
O*wy, a® /A*wy, Ge) Piles Vv sin 72% 
By? ay al Bx | | b2 erste Des vis vias P.n(4, hy) sinh nz c- 
By putting x = 7 uy: => Mand =) 0 ee ae + , the above expressions become 
l n—1 tanh" 
sblhe = F. p(s oO 
=| a eee) Die EO =) eae 4 
vO uN a oT Css) (1.3.5...) : 2 
St 
__ | wy, Me! CS Le) ee on (2 Meee eke eae 
ve ay? | a ecg Pa se, te 4 2 
aa ty te ears (e358) (315) Shee 
\ 
These expressions provide the following numerical results: 
ws. q a ats qe 
A = 0,116- " B= 0,038 sp: 


When the Poisson’s ratio is taken as y = 0.3, the bending moments can be evaluated as follows: 


M, =— D(A +B) =—0,0129qa@, My =—D(B+ yA) =— 0,00738 q a2, 


while the bending moments at the center point of the simply supported plate are as well known 
M, = My = 0,0479 qa?. 


_Accordingly, by superposing them, the bending moments M,, und M.,, at the center point of the plate 
in question are given as 


M, = (0,0479 — 0,0129) q a2 = 0,035 q a?, 
M, = (0,0479 — 0,00738) q a2 = 0,041 q a2. 


Finally it may be remarked that if a similar consideration about the antisymmetrical case of 
bending are made, a general case can be easily discussed by the method of superposition. 


9. Illustrative Example 2: Symmetrical bending of a simply supported 
plate clamped along the middle portion of four edges. If each two oppo- 
site edges are partly clamped in symmetrical arrangement (see Fig. 7), i. e. 


zy and moreover load distribution is symmetric, it can be written in (4.6) 
that 
ay a, | Ow Dee eo Je F=f = Fe 
Fig. 7. Plate with each opposite 


edges clamped symmetrically at In the above, 


: it has been taken that i = 1 since the number of the clam- 
the middle portions. 


ped portion on each edge is one. Then Eq. (6.7) becomes 


2S) s E, £650) +F.(0)} ®n(by h) ®,, (0,1, 


r? + __n2 


eee, a Fie Pants by) Pues) = 57 By By, )- 0) 
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Next, consider Eq. (6.8). It can be concluded that Px=0 =—Pxa Owing to the symmetrical 
distortion of the deflection surface and then B,, = — B,, by (6.1). And, we may consider n = odd 


because a symmetric function generally can be represented by an odd series. From Table 1, evidently 
faa) + f(a) = —[ful0) + fn(0)] - 
Similarly, from the relation aj = a—a,, it can be derived that m,s = odd. Hence Kq. (6.8) 


comes to coincide with Eq. (9.1). Same considerations as before prove the coincidence of Eq. (6.9) 
and Eq. (6.10), and give the following expression: 


DS) 8 FS [si(0) + FiO] Byala, by) Oya ( a fy) 


4 2 I (ar 

ie 2S, $ E, ey os | oe eis AS 1) ®,,,(a, h,) = >, as Parla, hy) 2 (9.2) 
s n r r2 + = n*| n 

b2 


From the conditions a; = a — a,, bj = b — b,, h, = a; — 4, l, = b, — b, the following expressions 
are obtained as in the previous example (see (8.3)): 


etic 
sin = 
a 
®,, (a, hy) = — 2 e » when m-—-n=even, 
m? — ane n 
ern, 
sin 
a 
®,,,(a, h,) = — 2 is » when s-+n=even, 
5? —— n? 
a2 
ony G; 
sin 
a h va 
®, (a, hy) = — 2 re , when s-+r= even. 
AS tees 
a2 


Then r also comes odd. Furthermore, the expressions for @,,,,(b, 1,), ®,,(b, 4), ®,,(b, 1,) ete. are 
made by replacing a,, a, h, in each of the above expressions by by, b, 1, respectively. Next referring 
to Table 1, it is known that 


where | 2 
0, =f, —tanhf, (8, tanhB,—1), B,=—5,-> 


Oya, tanhio(o, tankia, — Uy" Oe 


In the case of a uniformly distributed load, w’ is given by (8.6) and then A, in this case also 
is written in the same form as (8.8). In addition to this, considering the symmetrical properties 
of such a load, we can again use for w’ the expression which is obtained from (8.6) by the trans- 
positions (a, b) and (x,y). Hence the expression for B, is also obtained from (8.8) by the same 
transpositions, that is 


A, = 78 nt On » On = %, — tanh a,, (a, tanha, + 1), 
2 (9.4) 
B,=—22%9,, Ga =P, — tanh B, (f, tanh f, + 1). 
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From this, both the equations (9.1) and (9.2) are rewritten as follows: 


en bh Ss Ss P, (a, hy)] [r Pn -(, 14) 
Ls s E, ee P,,(b, L,) Drnlb, 1) =i we ae s By In ig ) b2 53 z 
n mu” a ( aE a”) 
a 


oki s nm a 
(13, 5.) (e3;5.5) (1,855.5) (U8:50.)' (8.5205) 


2 4b? re ; 
a i >; en Pra n(b, L,) ) (9.1) 


(1.3.53) 


h 3 80° [nD n(b, LY] [r Pn (ay) 
wy Fe DE Puls by) Pnalts hy) + oe DoE Dy (e+ ey 


(hagse) (1.3: 5) (352) (1.3, 5.) (1.3.5.0) 


2qa? , 
By oe Pal Baye C2) 


Especially, in the case when the plate is square and the clamped portions are equal to one another, 
no 
“Bt 
(9.1)’ and (9.2)’, E, = F, is derived. Thus both equations become identical and the conditional 
equation for this case is given as 


é 8 r 
>: s 1p 2 a ®,,,(a, h,) Pn n(a, h,) =I = ‘, aS. Guan ACH h,) ®,,,,(a, 4 


(1.3.54) (1.3.8...) (1.3.5.0) (1.3.5) 


Le a= ba, = bh, = l,, it becomes that «, = feria and. then:@,, = 0,,'0, == On Hence, from 


=e esa) (= 13 oe. (9.5) 


This represents an infinite set of linear equations with unknown s E,, the coefficients of which are 
symmetric functions about s and m so that the matrix of such coefficients is symmetric. Now, 
the double series which appeared in (9.5) can be transformed into a single series. So granting, 
the numerical computation becomes far easier. The second term in brackets on the left hand side 
of (9.5) can be summed up with respect to r (or n). Let us consider this simplification below. In 
the first place, observing that the following relation holds: 


ist 1 1 1 1 

2 2 2{- 2 ay ay OH ery 2 > 

(r2 4 0 (os ms (nt+ ie mt fe — im? Tore [n+ ae) (r? + n2y2 
i 1 i i 


\ 


the series about r is transformed as follows: 


sin) % 

F 

S +n ®,,,(a, hy) = — 2 >; ———* he 

(iis) Pine same) (m*— = r) 
a eee (Lara) a2 


a\2 

i yi ie nea SS r sin (2% 
\9 = — —_+>—— s1n 

a nt) - eae a r2 + n2 a 


2 we dy pole Cp Lee 2 3 
(n + he (153, Se) eu he” (4. 850.) 
GINE 
2 i) 
sob ait See ee (9.6) 
aeeicnne (r? + n2)2 
n+ 55m Gis. 


Next, the use of the following summation formulas of trigonometric series shall be made: 


4 conel(5—E) 
Sy . Tt 
Piape CUNT a She ge 7 O<é< iL) (9.7) 


r 4 
(1.3.5...) cos os 
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cosh 24(5—8) 


r : n 
(1.3.5...) cosh ye 
E : : foes A 
>; (CS es 61 E sinh 7 cosh a A(5—8) 
(1.3.5...) 8 A cosh? — 2 
2 
L . 1 wA 
—(z—8 sinh nA (gz —£) cosh | (Oreo aul ye. (9.9) 


Dee at 1 
By the substitution of a, = ee (a —h,) and the above formula (9.7), the first term in the large 
parenthesis on the right hand side of (9.6) is concluded to be 


maa [1 1 h, 
Dy = PR he Me > r sin (1 AE hy E Saat 
7.  r2— m2 * 5 r2 — s m2 2 A 4 aD s 
(1.3.5...) he (1.3.5...) he ee as hy 
cos 5 
Seechee (m = odd) . 
cos ce (F) 


And, the second term and the third term are respectively written by making use of (9.8), (9.9) as 
follows: 


cosh ers 
> r Be GB. BOE 2a 
: r2+ n2 Ga mead cosh 2 
GLEE) Toh 
r 5 PEE mv aah HOE paglin Wt so DVGE nat 
s = sav Mal eee 
> tn ine sinh 5 cosh a - sinh a cosh 5 |: 


nw 
(1.3.5...) 16 n cosh® —— 


If these summation formulas are used, every series on the right hand side of (9.6) may be rewritten 
in the form of a sum. Thus it is seen that the double series on the left hand side of (9.5) finally is 
transformed into a single series about n. Hence we obtain the following equation for actual com- 


putation: 


Nata ~ WGR ES 
, giv 2 : sin 
On a ox a 
Pas ig? a: It AE hea Ww tae Ki 2 oy 2 a Le 2 
Pi, 3.5...) Gis yee es ght cin, Pal (body soe ng 2 ie tgs 
, h 
: bs cos auth f Dosh e mK sink & 
2a nz 2a hy 2a 
x + i tanh = ES 
(m* coe n*) cosh > cosh —— cosh 5 
. nga, 
q a8 0, aa a A ] 3 c ) (9 10) 
x Qn Sel Oa eases a : 
— eh oe aU al Ave : : 
(ison tc ae gt 
where 
n nz (nz n9% 
On =—=-—t h 5 ("3 tanh 5 | 1), 


Now, suppose a case where every clamped portion occupies respectively the middle section of 
one-third of the edge, i. e. 
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Taking now the first twelve terms on the left hand side of Kq. (9.10) and calculating the coefficients 
with due consideration of the relation (8.9) in the previous example, we obtain a set of linear equa- 
tions shown in Table 8 corresponding. to m= 133.5, 2:25: 


Solving this system of equations, we find 


E, = — 0.263541 K , 
3 E, = — 0.354495 K , 
5 FE, = — 0.460605 K , 
TE, =— 0.565295 K , 
9K, = — 0.668501 K , 
11 Ey, = — 0.771339 K , 
13 E,, = — 0.873405 K , 
15 E,; = — 0.975489 K , 
17 E,, = — 1.07905 K, 
19 E,, = — 1.18236 K, 
2) Bo — 281532 Ke, 
23 E,, = — 1.38701 K 
a SEE where 
0.2, y-"eJa qa 
z0 z es 
05M, J (la 
ce y=Wlda From (7.5), the deflection function in this case 
ay | 7) is now expressed in the form 
1.0W, y=(lidla y=iza 
Wy=o.001g 22 iB ek HNL 
Fig. 8. Contour lines of deflection surface for the case when every in which (9.11) 
middle portion of edges is clamped. w"” = we, + we, + wy, + Wy, - 


Here, the deflection function of the simply supported plate w’ is given by (8.6) as in the previous 
example, and the remaining four terms are expressed by (7.1) to (7.4) respectively. That is 


l 


Me, ts, = — pa Dy Lhale) + fle) sin “22S” 5, o,,(6, 1), 
ES) (1.3:3:..) 
h = ; 
yt y= pas) Di enV! g0dl san Sang oe D(a | 
(1. 3:5...) (1.3.5...) 


Since E, = F., a = b, h, = 1, in the present case, we get 


wl Wah Wa eb Wy 
h, = ; “oe 
Set Dias a, sk, DE P,,,(a, hy) [Efal) + f,,(x)] sin ae + [gn(y) + 8n(y)] sin nae , 
(OPM YPGE ia 


where f,,(x) + f(x), Sn(Y¥) + 8n(y) are given in the item of Case (3), Table 1, and the computation | 
of ®,,(a, h,) requires the consideration of the relations (8.11). Thus, evaluating, by (9.11), the 
ratios of the deflections of some points to that of the center point as before, we have Table 9 and 
Fig. 8 in which the equi-deflection contours are shown. 
The moments at the clamped sections are given by (4.6). 

At the center of every section, i. e. 


such moments become 


My, = My = M, = MM, — 5) ee ee 
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The computation of the above alternating series by Euler’s method as before gives 
M,, = M,, = — 0.05931 q a? 


In the next place, let us consider the bending moment at the center point of the plate 


fe ee) >] , observing that the following holds: 
a 


2 tie 16 2 
n—l1 
2 
ow!” 1 aw!’ 2 hy 1) 
SS eRe see | Bars = sE Acie ae 
a Ox? a b ay? a b Daa 2 : 2, cosh 22 
sia Mele tiny Fe Dares O) (esr) (1. 3. 5...) 5 


qa 
xX @,,,(a, h,) = 0,0184 ae 
and taking the Poisson’s ratio y = 0.3, we can write the bending moment in the interior due to 
the clamping moments along the edge as 
M, = My =— D(A + 0,3 A) = — 0,0239 q a?. 


Therefore, the required moments M, 
of superposition as follows: 


M = M, = M, + MZ = (0,0479 — 0,0239) q a —= 0,024 q a2. 


» M, of this plate at the center are given by the method 


10. Hlustrative Example 3: Symmetrica bending of a simply supported plate clamped in the 
neighbourhood of every corner. If each two opposite edges are partly clamped in symmetrical 
arrangement (see Fig. 9), i.e. 


(10.1) 


Lely aah Inet 


hy =h, =h, =h,=h, 


and load distribution is also symmetric, the clamping moments must be 
distributed symmetrically on the two opposite edges. It can be, then, 


written in (4.6) that 
; tent ee oy 2%. 
s s,1 $,1 s (10.2) 


| Dp ed ree hE) teamed ey = 2g 


s 


Accordingly, such moments on the edge x = 0, for example, are expressed as: 


Fig. 9. Plate clamped near ever + SI - SH 
rene th Vi : E, sin 7 , Vie > E% sin = ; 
Ss s 
(1.2. 3...) (1. 2. 3...) 


The case considered now is such that the distribution of the clamping moments is also symmetric 


about the center point of one edge. From this, it can be taken that M,, = M,,, if y, =h —y;.- 
That is 


Ms, — Mz, = 3) E,sin 25%! S" ey sin S2™G—W) = Dues 1) Ei] sin $798 — 9, 


Hence we have 


and similarly 
Fe =(— 1-1 F% or F, = (— 1)*—' F, 
Next, owing to the symmetrical properties of the load distribution, slope displacements come to 


be distributed in a symmetrical manner. It can, therefore, be concluded that n = odd in (6.2) 


to (6.5). Further, since the relations; P=) = — Grea Pye = — Py=b must hold, it is easily 
seen that 


Bee oR, A, Ae (n = odd) (10.4) 
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Let us write down the actual expression by (6.7). Considering that the number of clamped sections 
is two to every edge, i.e.,a =f =y = 6 = 2, we have the following two equations: 


ie i 
FD) 8 Baa Dd} fil) Pyn(bs 4) On nlbs Ls) +2 S's Ex >) f5l0) Pen(bs bo) Pran(bs fi) 
ra = i = T if os 
does a 8 Ey >) frl0) Byu(b, h) Pn n(Bs 1) +2 >) SE, D7 FiO) Bulbs &) Pyun(bs U) 


a5 See sFiy oe SS 


n 


r ; 
b2 \2 ®, n(@, hy) Pn r(b, A) 
(" + — n' 
az 


4h, b 
tr Fas >, (en Mp Ponts a) Ppp (bs I) 
s n r r2 + =; n*) 
= 


Spe 
s n r (t+ a 


a 


* (a, hy) ©, ,(b, lj) 


| 


4h, b? = (—1)"1 nr es 2 
25 ar ai > 2 Fae as i ae eae ®, (a; hg) Drrld, A) 


aio 
= os B,, ®,,n(b, 1;) a ne ee ea (10.5) 


The other equations (6.8), (6.9), (6.10) may also be written in the same form. We now proceed to 
investigate the form of the expression ®. From (10.1) 


b, = b, =0, Rie pea b,=b, =b—l, ee fae 


and then the following expressions are obtained from (4.8): 


®, (b, ,) = ®, (6 4) = (—1)*- 


nal (10.6) 


and then @, (0, In) = (— 1)"** ®, (6, Li) - 


A similar expression for @,,,, is also obtained from the above by changing s for m. Moreover, 
@, (a, hy), ®,n(4 hy), Pyn(a, hy), P,n(4 hg) and so on are expressed by similar expressions as before, 
Bf Be Ty> los ty, 1,, Lin (10.6) are replaced by a, hy, hp, hy, hy, h respectively. Thus, by taking into account 
the foregoing relations (10.1) to (10.6), rewriting (10.5) gives 


ESPs E,W Ua) + F501] Paolo hy) Palo 
4h B? n—CKpyi—C rn 
ee DP sie en 


= >} By Paun(bs I) Gin anya Ml pee (10.5)’ 
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or 


21 0. 16h b2 [n ,n(a, hy)] [r Dm r(b, 1;)] 
ADM MOET a ct oS nee ) | 
) 


s r n 2 pa) 
CLO RHA eae ROR (1:2, 8..c)0b floes aise si) Ted goal 
= Or) (i 1285. ae ale (10.5)” 
(1.3.5...) 


And also, considering the relations (10.4), (10.6) and 
Ina) + fala) = — [f:(0) + fr (0)] = — 2 


Kq. (6.8) also is reduced to Eq. (10.5) in the same manner as Eq. (6.7). On the other hand, the 
same consideration as above gives the following equation with the coincidence of both the equations 


(6.9) and (6.10): 

2h 7 iG 161 a? aM Pn(b, 1,)] [r Gm (a, hj) 

ra taps a Ponlts by) Pala, hj) + ES S's E, SD) Det Dieter 
) 


eo, 2 2 
(1.3.5.0) cd 52) ( 1 5 oe, 


(E2is2) 6 ed (1.2.3... 


= >) An ®un(a, hj) (Ti 3 ea ef eal aoe (10.7) 


(ease) 


where @,, Q, are given by (9.3) in the previous example. Though Eqs. (10.5) and (10.7) alike re- 
present two cases when j = 1 or 2, these cases again come to coincide with each other by means 
of the relation (10.6). Conditional equations for this example, therefore, can be represented by 
the two systems of equations which are derived from (10.5) and (10.7) by putting hj =h,1l =|. 


Especially, in the case when the plate is aquare and the clamped portions are equal to one 
another, i.e. a = b, hi =1,, we have, referring to (9.3) 


ee hee Oe Pane aL er Rp eae Ls 
Ores Oh eae orn 2 ). 


If such a plate, moreover, is subjected to a uniformly distributed load, the intensity of which is q> 
the expressions already referred in (9.4) become 


9 3 
2qa 
eB a iB ape On» 
On = On = So tanh (2 tanh _ — i) ; 


Then the combination of Eqs. (10.5)” and (10.7) yields s = F, and accordingly both equations 
become one and the same. Thus the conditional equations for this case are reduced into one system 
of equations as follows: 


: Cn 8 
Br Fl Bs Peles) Pnltst) + SYD) aga we Balt by) Pe | 
( ) 


(1.2.3...) (3s) 1.8.5...) (e305. 


qa* in 
Ph, >) = P,»(a, hy) (m = 1, 2,3,...). (10.8) 


(U3n5:-.) 


Comparing this equation with Eq. (9.5), we find that the form of the equation is similar to the 
previous one and their coefficients again constituting a symmetric matrix of infinite elements. 
Here, the second term in the brackets on the left hand side of (10.8) can be reduced to a single 
series as in the previous example. Let us now consider this point. However, it will be mentioned 
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that this case differs from the previous example on the point where the index m here follows natural 
numbers, and about the form of @,,,(a,h,) and so on. First, we have as in the preceding case 


_ rah 
>) Gea Palas In) = (— 1)" Seas 
T CBS ae ee ) _ ae baa oe 
(1. 3.5...) (izs.5e)n. es (m* eye P| 
a\" 
orl 
a + as h a r rah r rah 
ie eet | ae eM 
h® (1.3.5...) h? (1.3.5...) 
a 2 
oft 
h \ r a pile 
h? (e355 20) 
and, by the formulas (9.7), (9.8), (9.9) 
mam ‘| 1 =| 
r rah cf ea ion a on a 
ee a 4 A eas SC AE iy Ue nay. 
B.D se he oh 
i hy 
cosh n a G —— 
Diseek a 4 — 2 ? 
(1.3:5...) zosh — 
ig ES ee me 1 Be = Gosh 1 h 
(r? + n?)? iam nn E vo ow na(a— a 
(1.3.5...) 8 n cosh*- 2 


1 beer 1 h nm 
—(3—G) sinh na (3 «)eosh ake 


Making use of the above expressions, every series on the right hand side of (10.9) can be summed 
up in closed forms. Therefore, referring to the expression of (10.9) thus obtained, the double series 
in the brackets of (10.8) is reduced into a single series, so that we obtain 


A sin? —— = pee 
___ 7)s—1 On 2 we | —— ae Le — 
aaa ee a a oe 
(1. 2:3...) (1.3.5...) (™ —an)(s Fite! (1.3.5:..)\° —3" UL +a" 
2 
n( cosh nz ¢ — +) coshna( 4 —*) 
x a 2 a spit pes 2 a 
oS = te a - | tanh —- = 
(m* a Hf n cosh ‘ cosh > 
a 
h _ nah 
1 h sinh nx( 5 = ete nsin~ g 
ees 12 = ; 
2 a i a - ‘ li 4 
oe (a apm ree | RLS aster 
h 
p See = 
eae See Be (m= 1, 2, 3,...) 
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In the coefficients of (—1)*"" s E, on the left hand side of the above equation, the series 


nah 


n sin 


> ee he 
(1:305:3) (m: Wea n ( Vig? n' 


again appears. This series, however, is of the same form as (10.9) and then can be summed up in 
closed form as before. The foregoing equation, eventually, is rewritten in the following form for 
actual computation: 


5 AND : 
sin? sin 


ah 
= On a a 
Dice! ” ie pe ey he h 
s n Ss —s. 


2 
Ph aa 7 a Pap Pari Aare) 
(1.2.3...) Cae) me) n (m a” (1.3.5...) (> oe )(m 3 a” 
: aes | 
an(*) cosh nar{ + — 4) 3 f costinx(p— 4) 
x + = tanh “ 
(m + => n*) cosh be d cosh ee 
2 2 
sinh nx(  — i 
I I 2 a Li (Bae ae 
x( =, cosh 4 m (s? + m?) h 
oe coe Came dale cont (FE 
h 2 a Che nh im h 2 a 
x (s? + m?) osh 2% ap 2 2 2h osh 2 
2h 2h 
amz(1 
1 h vies h z— ") } (— 1 It 
Hey ama 2 (s2 + m2 
cosh rh 
5 HEI | 
AR 6 sin 5 
rye ni], he : Crna V2 SS) oe vs (10.10) 
(1.3°5...) (m Ser Gaus 


| 


Let us now consider the case of a square plate when every clamped edge occupies a one-third 
section of the edge from the corner, i. e. 


h l 1 
a=b SS SS SS 
’ a b 3 
Putting h/a = x, we have 
. nah 
oy 1 sinnax Tt Tt 
a : 
lim 5 lim ——, = cos m7. = (— ])™—1 —_ 
Belt oy le ie 2m 2m 
m?——n gg RE) 
*% a2 aaa n2 


for the case when x = m/n, while m/n = h/a = 1/3; from this $i =p ce odd, and so (—1)™! 
= 1. Therefore, we find 


5 LUMA 
sin = 
a 4 h m 
he aoe when — == = 
WS ne a n 
a 
. nah (10.11) 
sin h | 
Te s 
a — = when —— 
sa Ue n2 2s? a nm 
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as Bee Sa as CT tg hen Eee ROD ON RS PRLS Dec efe|feal Qi con St non SS reo) (Gn Sac SO con 
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The deflection function w is expressed from (7.5) as follows: 


w=w +w’, (10.14) 


4 =O | wz) | (w,, + W;,) -! (w,, | Wy) - (wy, ++ Wy.) 


where w’ is given by (8.6) and every term of w” by (7.1) to (7.4), and then the arrangement with 


h l 1 : 
(10.1), (10.2), (10.6), a =b,h=1, = 3 ete. yields 


b 
sin” 
Lee oe SG) ea <a [ile + fu(e)) sin =" +. (gy) + Bal) sin | 
(is203.2) (nat Seymore 


In the above, the expressions f,(x) + f,(x) and SAY) + 8n(y) have already been given in Table 1. 
Thus we can evaluate deflection w by (10.14), using the values of (10.12) or (10.13). The actual 
evaluation in this manner gives the following values at the center point of the plate: 


Wy = 0.00196 ae s when the values of (10.12) are used, 
Wy = 0.00195 ee . . when the values of (10.13) are used, 


so that we find that both the above two cases give practically the same value. Table 11 gives the 
ratios of deflections to that of the center point and Fig. 10 represents the equi-deflection contours 
as in the foregoing examples. 


Table 11. Distribution of Deflection 


ts 1/12 2/12 3/12 4/12 5/12 6/12 
x/a 
1/12 0.0873 | 0.1452 0.2019 0.2232 
2/12 0.2432 0.3451 0.4230 0.4524 
3/12 0.4037 0.5403 0.6342 0.6682 
4/12 0.5403 0.7009 0.8075 0.8441 
5/12 0.6342 0.8075 0.9207 0.9594 
6/12 0.6682 0.8441 0.9594 | 1.0000 


The clamping moments are also given by (4.6). At the center point of every clamped section, i. e. 


Sines Pree, oe Rete eR OY oe Ye 1 
h, h hy ii, lb l i i 


2 2 2 


we have 


M,, = M,, = M,, — dy; = Te MG = My, = Ty By ee, ee 


Evaluating such a series by Euler’s method of summation, the following results are obtained: 
My, = 0.0288 qa, when the coefficients of (10.12) are used, 
My, = 0.0283 qa2, when the coefficients of (10.13) are used. 


: 1 
For the point i = 1 » the expression (4.6) may be arranged as 


Me 0707 eae 


where 
Ae (Eris) (Es te En) 4 (Eee 
BH he Ea Bi ee 


In this case, such 4, B can be evaluated by the Euler’s method. Using (10.12) gives: 


A= 06 (3 a B = 0.054 f SS 
4 
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and using (10.13) also gives the same values. Thus, we get 


2: 
M,, = — 0.049 (3 = = — 0.0149 qa?. 


3 
And for the OS arie = 4 > the following representation and the numerical result may be obtained 


M,, = 0.70711 A — B = — 0.0478 q a’. 


0 
) 
0.2Wo 
0.3Wg 
0.4Wy 
0.5 Wp} 
0.6W%p 
0.7% 
0.8Wp 
0.5Wo) 
LOW, 


y 
ga. 
Np=0.0095 “7 


Fig. 10. Contour lines of deflection surface for the case when every corner is clamped. 


1 ¥y 1 
Next, let us consider the bending moments at the center point of the plate(7 SS Z| : 
Using the values (10.12), it can be written that 
n—l1 
aA / = oe 
= we == oe = | ae 2 SS: Ss E (Ses 2 ®,,(a, hy) 
— 1 x? a b Ce pea ~ Daa or . = h e 
pe ae Gio (sit), Ges a) oN ~ 
i ie 
= 0.0141 D 


Taking the Poisson’s ratio vy = 0.3, the negative moments My, My at the center point caused by the 
clamping moment of the edges are given as 


M? = Mi —— D(C + 0.3 C) = — 0.0183 qa? 


Hence, the required moments at the center are obtained by superposing the free moments M,,, M, 
on the above reactant moments, i. e. 
M, = M, = (M; + Mz) = (0.0479 — 0.0183) q a? = 0.0296 q a?. 
On the other hand, if the values (10.13) are used, we obtain 
M, = M, = 0.0294 q a. 
The practical coincidence with the previous result is here recognized. 


28* 
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11. Experimental evidences. To verify the theoretical results obtained in the foregoing example, 
some experiments were conducted. Aluminium plates are used as specimens, which are held on 


TLOURSS 


slime 
ape edges 


SALLE LLP AN TESS EEIES TID OD ODILTS PETER 


LTE RSS WZ MX 
MIE 


Fig. 11. Detail of supporting device. 


However, 


a steel frame by putting them between sharp edges for simple 
supporting and by clamping with additional steel blocks for 
clamped section (see Fig. 11). As the loading apparatus, a 
lever is used for point loading and an air bag for uniform 
loading. The actual experimental results are developed in 
Table 12 in which, to show the accuracy of experiments, the 
results for a simply supported plate and a fixed plate are 
presented together with their well known theoretical values. 

12. Conclusion. A method of handling such a rectangular 
plate having special mixed boundary conditions was pre- 
sented and also its reliability was investigated. The labour 
required for computation amounts to a substantial quan- 
tity and thus an approximate method seems to be desirable. 


such a difficulty could be solved through the use of a high speed digital computer at 


present. Actually, the author has taken advantage of such a devise. 


In conclusion, the assistance of Mr. S. Hatano and Mr. 


H. Okamura on numerical computation 


and experiments is gratefully acknowledged. 


Table 12. Experimental data of the deflection and the strain at the center point of square plates 


Connie Deflection per Load 0.01 kg/em? Computed : Strain per Load 0.01 kg/cm? 
deflection Computed | Measured Diff. Moment Computed | Measured Diff, 
) 1* mm mm % q@ 10-4 10-4 % 
ee EE eee 
*1 *1 
Mx, = 0.0478 0.279 0.260 + 6.8 
0.0443 0.185 0.170 + 8.1 
My, = 0.0478 0.279 0.260 + 6.8 
ae *3 *1 
M, = 0.0350 0.203 0.204 — 1.0 
| 0.0323 0.146 0.143 + 2.1 
pate = | My = 0.0405 0.272 0.279 — 2.6 
az *2 - *2 
M, = 0.0240 0.157 0.172 — 9.5 
0.0155 0.0712 0.0750 — 5.3 
= My = 0.0240 | 0.157 | 0.172 | 9.5 
= *3 *3 i 
VE —10L0294 0.183 0.200 — 9.3 
a 0.0213 | 0.0954 | 0.0975 | 429 |—*— 
Ber My = 0.0294 0.183 0.200 — 9.3 
* 4 - : *4 
M;, = 0.0229 0.143 0.145 —l1. 
0.0138 | 0.0625 | 0.0665 6A = 
| My = 0.0229 0.143 0.145 —1.4 
1 The numbers *1, *2, *3, *4 denote the following kinds of specimens: 
*J, Aluminium plate, 300 x 300 x 3.00 mm, E = 720000 kg/cm?, 
*2, Aluminium plate, 300 x 300 x 2.97 mm, EF = 673000 kg/cm2, 
*3. Aluminium plate, 300 x 300 x 2.97 mm, FE = 690000 kg/em?, 
*4, Aluminium plate, 300 x 3002.96 mm, E— 690000 kg/cm?2, 
ae (Computed yalue— Measured value) 
ae Computed value aie uae 
* It is assumed that Poisson’s ratio y — 0.3. 


(Eingegangen am 17, Marz 1959.) 


Anschrift des Verfassers: Professor Dr. Muneaki Kurata, Osaka 


(Japan), Osaka City University. 
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Spannungen in anisotropen kreiszylindrischen Rohren* 


I. Mitteilung 
Von K.-H. Miiller 


1. Einleitung. Bei der Herstellung von Bauteilen kann man durch geeignete Hinzunahme von 
Anisotropien eine giinstigere Materialausniitzung und damit oft eine gewiinschte Gewichtsersparnis 
erreichen. Holz und verschiedene Kunststoffe sind durch ihre Faserstruktur dazu besonders pra- 
_destiniert; bei metallischen Baustoffen ist die Anisotropie im allgemeinen sehr gering, kann aber . 
durch eine entsprechende Walztextur verstarkt werden. 

Unsere Untersuchungen gelten der iiber den Querschnitt anisotropen Kreiszylinderschale (Kreis- 
rohr), die tiber den Umfang eine beliebig verteilte, aber langs der einzelnen Mantellinie konstante 
Belastung tragt. Hier liegt, abgesehen von der Randstérzone, ein ebener Spannungszustand in 
einem Ringbereich vor. Wir wollen dabei einen etwa schon im Kreisring vorhandenen Higen- 
spannungszustand unbeachtet lassen. Die Gleichgewichtsbedingungen befriedigen wir identisch, 
wie iiblich im zweidimensionalen Falle, durch Einfithren der Airyschen Spannungsfunktion U(x, ¥). 


2. Integration der Differentialgleichung fiir die Spannungsfunktien. Die Hinzunahme des Hooke- 
schen Gesetzes fiir anisotrope Kérper 


Ex = $11 Ox + S12 Fy 1 S16 Txy > 
Ey = S12 Oy 1 S22 Fy 1 S26 Txy » (1) 
Vey = $16 Fx 1 $26 Fy 1 S66 Tay 
fihrt tiber die Vertraglichkeitsbedingung zur Differentialgleichung fiir die Spannungsfunktion: 


aU aU aU thal OF aU 
S99 pa 2 $28 Fd ay T (2 Siz + S66) G2 Gy? 2 $16 de ays + SU ay 


=0. (2) 


Die charakteristische Gleichung von (2) besitzt, wie Lechnitzki gezeigt hat, stets vier komplexe 
Wurzeln: t, = ty, ty = t,, wobei der Querstrich die Konjugation anzeigt. 
Die allgemeinste Form der Spannungsfunktion finden wir dadurch, da® wir (2) mit 
See a ep or ey (3) 
transformieren, zu 
4 
U= a ight ey) 

wie man leicht durch Einsetzen in (2) verifiziert. Da U reell sein mu, erhalten wir daraus 

U(x, y) = 2 Re [p&) + v)] » | (4) 
wobei 7(z); P(e) willkiitliche Funktionen ihrer Argumente sind, die gema den am Korper an- 


greifenden Randkraften spater festgelegt werden. 


Setzen wir zur Abkiirzung noch 
/ , d 2, 
Pp =— dpa) A = (zo) 
dz, dz, 


Pi = Su — Sieh + S22 Pe = $11 t5 — $16 bt 1 S12 > 
BaP tt (a en 2 
t 1 = S1at — Soe 41 + S2a> $2 Io = S12! — See ty Pr Spe > 
dann werden die Spannungen und Verschiebungen bei bekanntem pz) und *p(Z2) festgelegt durch: 


o, =2RKe [Bo +h", 


o, =2Kelp’ +y'], (5) 
Tey =—2Re lag” +ey'l- 

u=2Ke[pig’ +py)> 6) 

vo =2Relay’ + qo" | . 


* Herrn Prof. Dr. H. J. Seemann, Saarbriicken, zum 60. Geburtstag. 
1 Siehe G. N. Sawin, Spannungserhéhung am Rande von Léchern, S. 34, Berlin 1956, 
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Die am Rand unseres Bereiches wirkenden auReren Lasten erfassen wir mit 


4h, 
h 
hal Ya [Yee : (7) 
l 
I, 


1 
Xa [X,b— 7 
I, 


dy Ls 0 


Ks empfiehlt sich beim Kreisring, die radialen und azimutalen Spannungs- und Verschiebungs- 
komponenten zu benutzen. Wir erhalten sie aus den Invarianten zu 
Op +0, = Oy Oe 
Og — 6, 2htg — Oy Oped eats (8) 
v, t+ivg =(u +ivjye®, | 
Um den Spannungs- und Verschiebungszustand im Ring bei vorgegebenen Randlasten (o, — i Tr) Rand 
aufzufinden, entwickeln wir diese in Fourierreihen 


++ 00 ; + 0 
(([—U7,s),2n = DA, (6; = b3 4) io Se ee (9) 
Die Relationen (5) und (8) fahren aber zur Gleichung 
2 (6, —i t,o) = [C1 + 4) — (1-44)? 21] yo” + [(L + 8) — (1 Fit)? 249] ge” 
+ [Ql + &) —(L— itp ei] yp” + [(L + 8) —(1 + it) 2] y”, (10) 


in der also nach (3) (9) z, 1,2, als Argumente vorkommen. 
Um hier ein einheitliches Argument zu erreichen, veranschaulichen wir uns, da die Gleichungen (3) 


unseren Kreisring { der z-Ebene (s = « +iy =re’®) affin auf den elliptischen Ring RR. der 
. 2 


z ie ey av : ie ers 

komplexen | ‘bene | i Lee projizieren. Die zugehérigen Transformationsformeln 
% ig = X_ Tt 
2 2 2 2 


1l—i aoa = 1 —it I Mi 
= sts eS, So se | ae (11) 
sind aber nicht analytisch, da z darin auftritt. Wir normieren so, daB die Zylindermittelflache 
den Radius |z| = 1 hat. Fiir die Punkte jz] = 1 ist = 1/z; aus (11) folgt dann 


aT 


_ bev, Se aie 6 Sts lit, 
eee th wea a Sey ca Manoe a2) 
Bei einem Kérper von (beziiglich des Durchmessers) sehr geringer Wandstarke, einer Schale 
also, diirfen wir diese Abbildungsfunktion fiir den gesamten Ringbereich benutzen. In Gleichung (10) 
1a8t sich dadurch die einheitliche Variable z verwenden. 
Wir setzen nun: 


2 


ee) = Saye". w"(ale)) = 5 by 2" (14) 


in (10) ein und erhalten durch Koeffizientenvergleich entsprechender z-Potenzen zwei Rekursions- 
formeln fiir die noch zu bestimmenden Koeffizienten a, b,. Zwei weitere Rekursionsformeln folgen, 
wenn wir n durch — (n — 2) ersetzen und anschlieBend konjugieren. Somit stehen zum Errechnen 
der Koeffizienten die vier linearen Gleichungen zur Verfiigung 


R*(L+ a, +R" (0+ a) o_,— R21 ina, , 
= RL) (Ea pi eh ee 
oe" (1 ei 1) a, =e On (1 a tt) a, aa ee (1 ox t,)° a,» 
eras me (1 aig t,)? 4_(—2) Tipe rac [b,] =2 Bio. ? 


15 
—R(1+ 4) 0,—R" (ita, + R14 a)a,_, 2 
ot) (ei) A(n_z +++ [b,] = 2 Ane. OS), 
a (lat by aan — Oe (ith Guest eee Fy aes 
+ o*-9) (1 + 2) GS Say tee al ee Bs, on 


Da sich einer vorgegebenen Belastung nur eine bis auf einen iiberall konstanten Grundspannungs- 
zustand eindeutige Spannungsverteilung zuordnen 1aBt, kénnen wir a, by willkiirlich wahlen: 


a — bi Or (16) 
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Das Einbeziehen der Randlasten nach (7) und die Eindeutigkeit der Verschiebungskomponenten 
bringen 
; ——— ; ——— + Xi 
(@ +in)e, Oia, +0 tigt,—O—im)b, =A 


(py tig) @-1— (Pi —* hs) a_, + (pz + 7 qe) b_, —(p2—igq.) 61 =9, 
—(1—it)a_,+ (1 +i) a_j—(1—t4) Nie Gey Gb Bel A Wet aa en 


— (py, —tq;) 41 + (Pr + i qy) ay —(P2—# 4p) 6-1 + (P2 + i qo) by = 0. 
Stehen alle am einzelnen Rande wirkenden Krafte im Gleichgewicht, dann ist a_, = b. =0, 


bo 


Die Rekursionen (15) bilden, wenn man die Koeffizienten . b 
—1 9-4 


heranzieht, jeweils ein System 


von vier linearen Gleichungen fiir a, by 4_2b_» sw. Da die Koeffizientendeterminante von (15) 
@ 9; Gg 9-8 

hierbei aber proportional der Wandstarke e ist, und wir mit einem e < 1 arbeiten, wird man zweck- 

maBig das System etwas umformen mit 


R" =1+m-—, om =1—m> (18) 
zu ee 
@ 4 tye, $0 +8) e5, 0 —i4)8o, 2 — CO Piao gg eo = G), 
(+) o,—O Fa) a, — 2 (1 — it) ag +2 OF EP yt = OP, 
(1+ 47 o, +014, +8) 4,.—C FB) ogy $= GD), oa 
(4 ino, —GC—iafa, ——" 4 Bat OF Dogg t=O, 


wobei 


CoA. aa B,, ai (A,—B,)n-> , 


Cet aD a Ce 
pas 

2 

i « ES = e 
C® = — (Ans) + Bn») + (An—2) — Br) @—2) = > 

Be ee i B oie 
Cc i —(n—2) ew =2) { (A_(1—2) aE B_(n—2)) n 
12 


2 
ist. Damit ist die Koeffizientendeterminante aus der Abhangigkeit von e befreit. 
Bei einer p-mal stetig differenzierbaren Randlast sind die Fourier-Koeffizienten A,,; B,, vou 


der Ordnung 0 (sera): Man kann also hier mit der Substitution: 
m2+ & 


On 
O, = Tap? bere 

die Konvergenz der Reihen 7 und y giinstig beeinflussen. 
Unsere elastizitatstheoretische Aufgabe ist damit auf die algebraische, das fortgesetzte Lésen von 
vier linearen Gleichungen zuriickgefiihrt. 

Aus den so bestimmten Koeffizienten erhalt man ¢(z,(z)) und w(z(z)) und damit die Spannungen 
und Verschiebungen dadurch, daB man die Abbildungen z, = 2,(z) und z, = z,(2) riickgangig macht. 

Uber (12) finden wir 2(z,) und 2(z,) und erhalten so 


sees Si Ee 
oa) = Dio, ATED), 


i 1—ith (20) 
* oy ea eee 
va) = 3,4, (2+ ae ee) 

4 


3. Einzelkrafthelastung. Bei stiickweise stetiger Randlast bedarf diese Vorgehensweise keiner 
weiteren Exlauterungen; bei Einzelkraften wird man zweckmafigerweise den folgenden Weg 


einschlagen. 
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Uber den einen Rand unseres Ringes seien N Einzelkrafte P, verteilt, die mit der Last des 
anderen Randes im Gleichgewicht stehen. Wir ersetzen zunachst jede von ihnen durch eine iiber 
den ganzen Umfang stetige und einmal stetig differenzierbare Funktion, z. B. 


0 fiir 0<0Sc¢,—aJ/a, 
FO) = 4 k, [1 + cos(ad + .¢,)] fir ¢,—ajaS9Se,+a)/a, (21) 
0 fir c,t+afaxsd0S2n. 


Mit einer solchen stetigen Randlast (« wird positiv ganz und so groB gewahlt, daB sich die einzelnen 
f,(9) gegenseitig nicht tiberschneiden) errechnen wir nun die zugehérigen Spannungen und Ver- 
schiebungen. Danach fihrt ein Grenziibergang 


la ep (22) 
ky, a—>0o ee 


zur gewiinschten Spannungs- und Verschiebungsverteilung bei Einzelkraftbelastung. 
Hine solche Darstellung ist aus zweierlei Griinden der oft benutzten Rechtecklast: 


| 0 
f(P) = 4 & (23) 


vorzuziehen: 
1. Ihre Fourier-Reihe ist im Gegensatz zu der von (23) fiir jedes endliche a absolut konvergent. 


2. Da ja tatsachlich niemals ein linienférmiger Lastangriff (Schneidenlast) auftritt — Her- 
stellung und elastische Deformation von Schneide und Werkstiick sorgen dafiir — vermag sie eine 
konzentrierte Belastung sicherlich wirklichkeitsnaher zu erfassen als eine Rechtecklast. 


4. Bemerkung. Die Auswertung unserer Ergebnisse fiir technisch-wichtiges Rohrmaterial kann 
erst in einer spateren IT. Mitteilung erfolgen, da die Messungen zur Ermittlung der Anisotropie- 
konstanten von Metallen und Kunststoffen noch nicht abgeschlossen sind. Die in einem weiteren 
Bericht folgenden Untersuchungen erfassen dann auch im Rohr vorhandene Eigenspannungen, die 
wir in der vorliegenden Arbeit unbeachtet lieBen. 


5. Zusammenfassung. Unter Hinzunahme komplexer analytischer Funktionen laBt sich der 
Spannungs- und Verschiebungszustand eines iiber den Querschnitt beliebig-anisotropen Kreis- 
zylinderrohres unter willkiirlichen zylindrischen Randlasten in Form unendlicher Reihen be- 
stimmen, deren Koeffizienten rekursiv aus einem System von vier linearen Gleichungen errechnet 
werden. Eine geeignete Korrespondenz zwischen Anisotropierichtung und Last ermoglichen eine 
optimale Materialausniitzung. 


(Eingegangen am 9. April 1959.) 


Anschrift des Verfassers; Dr. Karl-Heinz Miller, Institut fiir Technische Mechanik, 
Universitit des Saarlandes, Saarbriicken. 
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Eine Verallgemeinerung des Prinzips vom Minimum 
der potentiellen Energie elastischer Kérper 


Von D. Riidiger 


1. Einleitung. Die aus dem Prinzip vom Minimum der potentiellen Energie abgeleiteten direkten 
Methoden zur numerischen Darstellung der Lisung eines vorgegebenen Randwertproblems der 
Elastizitatstheorie haben fiir die praktische Durchfiihrung den Nachteil, daB Funktionen zu finden 
sind, die auf bestimmten Teiloberflachen des Kérpers vorgeschriebene Werte annehmen miissen. 
Durch eine Modifikation des Prinzips vom Minimum der potentiéllen Energie gelingt es, die Rand- 
bedingungen des Variationsproblems fortfallen zu lassen. Dadurch wird es méglich, bei jeder vor- 
gegebenen Randwertaufgabe der Elastizitatstheorie mit — hinsichtlich der Randbedingungen — 
willkiirlich gewahiten Funktionen die Ritzsche und mit beliebigen Partikularlésungen der Grund- 
gleichungen die Treffizsche Methode anzuwenden. 

Zur Ableitung des Variationsprinzips werden die fiinfzehn Grundgleichungen der Elastizitats- 
theorie kleiner Verschiebungen in allgemeinen Koordinaten 


o*4|, + p*=0, | (1,) 
Dim = = (v4 + we), (i,) 
o’# =26 pee + p58 D;) (1) 


verwendet. Hierin sind o?“, D’“ und g?“ die kontravarianten Koordinaten des Spannungs-, Defor- 
mations- und MaBtensors. Der Verschiebungsvektor ist durch v*, der Vektor der Volumenkraft 
durch p* gekennzeichnet. Als elastische Konstanten werden der Schubelastizitatsmodul G und die 
Querdehnungszahl y benutzt. Sind n, die kovarianten Koordinaten des Normalenvektors der Kérper- 
oberflache, dann ist ihr Spannungsvektor durch 


ta ote ny (2) 
festgelegt. 


2. Das Variationsprinzip. Ein elastischer Kérper mit der Oberflache F unterliegt der Volumenkraft 
ip = p Ae = ye) und ist beliebigen Randbedingungen De Dabei sollen auf den Teil- 


flachen é, G, G die Flachenkrafte #, 2, es und ay gen, Teilflachen i 6, C die Verschiebungen v1, v”, v® 


vorgegeben sein. Es ist dann F = C + G, Pi C = ¢ und F = G + C. Zur Formulierung des Vari- 
ationsprinzips wird die Formanderungsarbeit —d// des im Gleichgewicht befindlichen elastischen 
Korpers bei dem Verschiebungszustand dv* gebildet. Sie ist mit dem Volumenelement dV fir einen 


Korper vom Volumen R 
— II = fff (o*|, + p") ov, dV =0. (31) 
R 


Fiihrt man in (3,) die Identitat 
o|, dv, = (o*# dv,)|, — 07" OD, , 
ein, dann ergibt sich 


—dlT = — fff o#dD,,dV + fff (0% dv,)|, dV + fff pou,dV =0. (3.) 
R R R 
Fiir den Integranden des ersten Volumenintegrals kann mit dem Hookeschen Gesetz (1g) 


o" 6D,, =2 G(Di SD DAD, _ 


eens yDiDh| 


bo| eH Be 


5 (0% Dj,) => 6 (0% m4, 
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(074 v,)|,, - p? UV, > 
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geschrieben werden. Beachtet man noch die wegen (1,) giiltige Beziehung 


42,2, 
oy, |, = 
1 
+0 (07 v,)|,, + p* v4] 


so folgt schlieBlich 
of OD; — 
Fir die potentielle Energie 6/7 erhalt man damit 
a= s[[ [eomephar +s {| | [oe Ie 
R k 
— | [ ferme nav—|f{ hatha (4,) 
R 
Die Anwendung des Gaufschen Integralsatzes 
PPC hav =f 0) naa 
auf das erste und dritte Volumenintegral in (4,) ergibt 
él = 5 [ focrenmart ys | { [aw Auie= 
F R 
— f fomoo.mar— | | [ 1$y,dV =0, (4,) 
R 
und nach Einfiihrung des Spannungsvektors (2) wird 
1 1 
=| foenats || fe p'v,) dV — 
F R 
| — | [ronan [ | [pena =o. (4,) 
R 
Leas 1 a a} 
Ersetzt man die Oberflachenintegrale in (4,) durch die Gebietsintegrale iiber G,G,G und G:G-e 
so folgt 
ae ee sce dfs eg 
sag LF fae) rit foe if L fay fe 
: : é 
aes ges 
_| [emma J feonar— f fesnars 
2 a 
G co 
(44) 


|! 


[fy Leo ay ae 


+a | [ [eed 


—t dv, = — o(t' v,) + df», , 
2 


—t? dv, = — 0(t? v,.) + Of v,, 
— 0(# v3) + dt? v, 


—t# dv, 


Fiihrt man in (4,) die Identitaten 
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ein, dann erhalt man zuniachst 


ee af+ + ie oe Byne ee 
wee Lf fae ) df +> / (2) df +5 [foe oe 
| dv, df — oF 2 dv, df — [ [ ®ssar— 


[rma 


3 | Mite.) dp II? eeu eee df + 


Ol|- 
Ae 


Ale 


D9, 


cel 
a [ fornars [Joa [Jona 


ae 2 | | [oer eerie Lip d= 0 


und nach Zusammenfassung 


i nea: ee Hears 
E [Jeena [Jarwa—t |frwe— 
= [fina [fos [rma 
+ [formar Jarwars foenars 

A ffamser[j fomer—v. 


al 
ales, 


alle 
ke 


Q|# 


eae 


(5) 


Die Beziehung (5) ist eine Modifikation der Energiehauptformel, aus der sich sofort das Prinzip 
vom Minimum der potentiellen Energie elastischer Kérper ableiten laBt. Da ERAS auf den Teil- 


Sar Ie 


oberflichen G,G,G die Flichenkrafte 2, 2,8, auf den peusherichen 6.6, é, é die Verschiebungen 


>}, v2, 0? und fiir den ganzen Kérper die Volumenkraft p p’ fest vorgegeben sind, gilt 


a= 8/5 [[e—2Ana+z [[@—2Ama+z | [C—2 ad — 


Bel Bs 
G G 


é 
a oe eae | ee aah 


a» 
Alle 


nee mG —2 p’) nav] =0. 
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Hieraus ergibt sich, daB unter allen méglichen Verschiebungen eines elastischen Kérpers, die auf 
der gesamten Oberflache beliebige Werte annehmen kénnen, diejenigen eintreten, 


die den Ausdruck 


nhs [Jeo See [Je »,6f — 


3 


1 ffee Sn hee a (eaeeo) ‘eee Pe Reon. 
— || [fz rn (6) 


zu einem Extremum werden lassen. Dabei hat das Extremalprinzip (6) gegeniiber der klassischen 
Formulierung! den Vorteil, da an die infinitesimal benachbarten Systeme von Verschiebungen 
hinsichtlich der Randbedingungen keine Forderungen gestellt zu werden brauchen. Man kann also 
zur numerischen Darstellung der Lésung v4 = v4 (x1, x”, x°) eines beliebig vorgegebenen Randwert- 
problems der Elastizitatstheorie im Sinne von W. Ritz einen Ansatz 


Oe 
Dice 


B= Seg (1) 
v=] 
mit willkiirlichen Funktionen gy machen. Hierbei sind die HUE aee g so allgemein auszuwahlen, 
daB bei hinreichend groBem n und passend gewahlten Koeffizienten '¢ ’ die Lésung der vorgegebenen 
Randwertaufgabe mit beliebiger Genauigkeit dargestellt werden kann. 
Das System von Verschiebungen (7,) liefert fiir den Deformationstensor, den Spannungstensor 
und den Spannungsvektor 


Bund SUlgh a vp 


vol 


Als Volumenkraft erhalt man 


piel) Git —2G (De ae ae git D:), 1 = Ge Ny. (7,) 


p! ca ae ot, D (73) 


Die Koeffizienten "t’ (vy =1,2,...,n; A =1,2,3) werden nunmehr so bestimmt, daB die mit (7) 
gebildete potentielle Energie ; 


Has{/ b— 2/8) baf+— >| { b—22)d, df-+> [J (® —28)b,af— 


1 2 
cel 


is soy ses ee Sea 4 [6 $22 oa ldpee 
ad mK Pp 4__2p \y, dV (8) 


dem Extremum miglichst nahe kommt. Die Bedingungsgleichungen fiir das Eintreten des Extrem- 
wertes lauten dann 


Q| 


Qo 


n 


all 
oan ee ae ee ee og} (9) 
c 
und haben die Galerkinsche Form. 
3. Das Verfahren von Trefftz. Beim Treffizschen Verfahren wird der Ansatz 


v= > cy (10) 


2 al. Zs B. C. B, Biezeno u. R.Grammel, Technische Dynamik, 2. Aufl., S. 80. Berlin/Gottingen/Heidel- 
erg 
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verwendet, bei dem die Funktionen q Partikularlésungen der homogenen elastischen Grundgleichungen 
sind und entsprechend dem Variationsprinzip (6) auf der Oberflache beliebige Werte annehmen 
kénnen. Bildet man die dem Verschiebungszustand (10) entsprechenden Elemente des Deformations- 


n n Lis : : 
tensors D’“, des Spannungstensors o*“ und des Spannungsvektors t’, so ergeben sich mit 


dtu = + (9 


ue + oe 4), 
aA r, r la A ya 
gw = 2G (die + 7a, gle da) 
und t” = s’“n, die Ausdriicke 


i Beer Bhs, Bh a HN mB 9, 
Dis Dy CO! Ge Sess Coe re tee (LS) 


ell y=l1 v 


. : v . “A ° 
Da die Funktionen gy’ voraussetzungsgema8 Partikularlésungen der homogenen elastischen Grund- 
gleichungen sind, ist 


Pal, 20. (11,) 


Die n linearen Gleichungen fiir die Konstanten lauten dann 


OM sels Drea 
oc 
und man erhalt mit (141) 
Mr =0= SSS ( by +g, —209,) df + [S(@e, +2 G,—28 pe) of + 
6c a 


rel 


S(#i,+89,—20 ps) op — 5 SS (28, + 8G — 283) 


3 


a| 
aH 


—fffpgaVv (is 1, 2,)5 5 «5 Mt). (12) 
R 
Bei Anwendung des Bettischen Theorems* 


Ue eee Peed Dl haat lene tone 


Gi 


= ppBaars [faa [Phu sbi fe Ji 
é Ci e © 


Q|- 


Me 


auf (12) kénnen die n linearen Gleichungen fiir die Konstanten ¢ (v = 1,2,...,n) somit in der Form 


nee (is — ®) G5 of — 


[fsa fr _ safle Bat 


i 
ine p,dV = 0 (vy =1,2,...,n) (13) 


angegeben werden. LaSt man in (10) die Anzahl der Glieder iiber alle Grenzen wachsen, dann liefert 
(13) die Koeffizienten des gegen die exakte Lésung konvergierenden Ansatzes. 


1 Vgl. z. B. C. B. Biezeno u. R. Grammel, Technische Dynamik, S. 87. 
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4. Die de St. Venantsche Torsion als Beispiel. Auf das Beispiel der de St. Venantschen Torsion soll 
das verallgemeinerte Prinzip vom Minimum der potentiellen Energie (6) angewendet und fiir 
einen geeigneten Ansatz (10) die Herleitung der zugehorigen Trefftzschen Gleichungen gezeigt 
werden. Dazu wird entsprechend Abb. | ein prismatischer Kérper in allgemeinen Koordinaten x*, 


Abb. 1. Geometrie. 


x3 (x, 8, y = 1,2), die mit den kartesischen Koordinaten x‘, 2 (i,j,k = ie 2) durch a! = ? (<1, x7); 
x3 — x3 verkniipft sind, beschrieben. Die Tangenten- und Gradientenvektoren ergeben sich zu 


k = 
k _ O% pelea tt te 
Ciena) ¢=c,=—90, ejta=T 
und 
ox” 
er Briel oe Yaak 
kT axk? Gq, = ce = 0, eal 


Die Randbedingungen der freien Torsion lauten in kartesischen Koordinaten! auf der Mantelflache M 
t—F —0 
und auf den Stirnflachen F’ und F”’ 


= Tay iy pus. 
Vj = — HE, XI x, ia—0e 


wobei ¢;; den ebenen ¢-Tensor und « die Verwindung darstellen. In allgemeinen Koordinaten wird 
analog auf der Mantelflache 
tb 0) (14,) 


und mit y® = cf «* auf den Stirnflachen 


V, =—KEgy? Be, P=0 (14,) 


pare selieben, AuBerdem ist die Volumenkraft p+ = 0. Fiir die in (6) auftretenden Teiloberflachen 


oe Be 
G, G, G und G, G, G gilt also 


3 


i oe a 
C= M C= GM orn 


1 Vel. z. B. C. B. Biezeno u. R. Grammel, Technische Dynamik, S. 120. 
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so daB die potentielle Energie fiir die de St. Venantsche Torsion in der Form 


M=z{ [2nd +z | [ee ad+g | | @—27) nat 


role 


+z] [@—2 edz | fe —26)4— 


F’ 


—z |] [{@—2Pyaa (15) 


angegeben werden kann. Fiihrt man in (15) die Randbedingungen (14) ein, dann erhalt man 


H=7[[edtz[[eadts | fenud+s | [ea 
M M 1a BY 


1 od , il fod wt 
Se (M5 +20 egy" 2) af —y | | ley + 2a eapy? a) df + 
BY Fv 


+z f{ {rear 
R 


Wird nun entsprechend (10) der fiir die reine Torsion zweckmaBige Ansatz 


ep (16) 


1 


[| * (v, —23,) df — 


na BB ee n 
V_ = —AEqpy” x, Vg = 49 =O 


IMs 


v 


verwendet, wobei p = p (x1, x2) Potentialfunktionen sein mégen, dann ergeben sich fiir die mit (16) 
gebildeten Elemente des Deformations- und Spannungstensors 


D6 = D® =0, Des (lees + ql), 
=| Ye g\*) (17) 


o*8 — 633 —(), G3 = Ga (—e*? yz + GI") 
und fiir den Spannungsvektor 


pan ngs 
(= 07" Ns. PEN ON pee ; (17,) 


. . v . . . 
_ Wegen der vorausgesetzten Potentialfunktionen p = p (x1, x2) ist p? = 0. Die in (16) auftretenden 
Konstanten ¢ (vy = 1, 2,...,n) werden so bestimmt, da die potentielle Energie 


n cua nn il non 1 rt 
Haz [fears | [Pharrs [Pears | fr ia— 
M M Me 


da 


ein , 1 OTRO uw 
=F [ [PC t2aeore)d—z | [PC +2ae%e ) df (18) 
F’ Les 


dem Extremwert méglichst nahe kommt. Mit (16), (17), F =i" = F “und | = 23” — x* erhalt man 


n n n n l n nt 
T= 5 [ a? mapas dee —F | [Bean x P= FZ pea (—erhyp + Hr) ma bdo — 
C1 F c 


l n 
—F | [65 etye +o) ano al- 
F 


Da é,5 y?|* = €48 g* verschwindet, folgt daraus 


HZ = pea (—e? y5 + H) mbas+ [| Carpytdf— | | Ca (eaoy" A) |= df. 
Cc " im 
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Die Anwendung des Gaufschen Integralsatzes auf das letate Integral ergibt schlieBlich 


= P60 (—e" yp +H) nepads+ 6a | [ ypxf— PCa tans? m pas 
Cc F Cc 


=~ 26agesay m pds + CaPsirnapds +a | | rex? df 
c F 


Cc 
n 


‘ A : oll 
Die Bedingungsgleichungen — = 0 lauten nun 
6c 


— 26a eggs? m g ds + 6a (Gl nef + $l m9) ds =0 aa bepp'e sears 1) 
é 


Cc 


und nehmen mit dem Bettischen Theorem fiir die reine Torsion 


Polar Gas — DP lan pas 
C Cc 


CaP Bla—eapy4) mi Gds=0  =1,2-..40) (19, 


die endgiiltige Form 


oder 


fd? nggods =0 G2 eee) (19,) 
an. 


5. Zusammenfassung. Es wird ein Variationsprinzip der Elastizitatstheorie aufgestellt, das als 
Sonderfalle das Prinzip vom Minimum der potentiellen Energie und das Castiglianosche Prinzip ent- 
halt. Da Randbedingungen des Variationsproblems nicht auftreten, kénnen die bei den Verfahren 
von Ritz und Treffiz verwendeten Funktionen auf der Oberflache beliebige Werte annehmen. 


(Eingegangen am 16. Mai 1959.) 
Anschrift des Verfassers: Professor Dr.-Ing. D. Riidiger, Freiberg (Sachsen), StraBe der Einheit 12. 
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Der erste Band iiber ,,Grundlagen und Anwendungen der Informationstheorie“ gibt 
einen Uberblick iiber die verschiedenen Zweige der informationstheoretischen For- 
schung im besonderen Hinblick auf die menschlichen Kommunikationspartner. Er 
beginnt mit der mathematischen Beschreibung von Signalen und linearen Ubertragungs- 
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lichen Sinnesorgane als Informationsempfanger leitet tiber zum Problem » Signal und 


Zeichen“ und damit zur sprachlichen Kommunikation in ibren verschiedenen Aspekten. 
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